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ABSTRACT 

Charcoal, charcoal residue, potting soil, aluminum foil, 
bismuth germanate, and petroleum samples have been 
investigated using instrumental photon activation analysis 
(i.e., no radiochemistry). The major and minor elements 
routinely observed by this nondestructive method were: C, 
Cl, Ca, Fe, Mg, Si, and K. A comprehensive review of the 
principles of IPAA was also included in the study. The 
principles were applied to a theoretical analysis of an oil 
sample in which the trace element concentrations were known. 
It was concluded that IPAA is a highly sensitive technique 
which could be used to fingerprint oils. 
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I . INTRODUCTION 



Activation analysis is a widely used analytical procedure 
for the study of elements in a matrix. In most studies, 
thermal neutrons have been used because of the general 
availability of nuclear reactors and the large cross sections 
of many nuclides for thermal activation analysis. However, 
not all elements can be conveniently determined by neutron 
activation analysis. For example, carbon, nitrogen, oxygen, 
iron and lead are not highly activated by thermal neutrons. 
In other cases the (n,tf) reaction may result in a product 
which is stable. For these situations photon activation 
analysis may prove useful. In fact, photon activation offers 
several advantages as a complement to neutron activation, 
(a) As alluded to above the (Y f n) , ( 5f,p) , et cetera, 
reactions of target nuclides often lead to products other 
than those resulting from neutron activation. Thus there are 
more possibilities of forming products with half-lives and 
gamma ray energies convenient for analysis. (b) For a 
matrix which contains elements with very large thermal 
neutron cross sections, flux pertubations can cause errors in 
the procedure. This problem, known as self-shielding, is not 
as serious with photons as photons typically have longer 
ranges in samples than neutrons. (c) Finally, since 
photonuclear cross sections are a function of photon energy 
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it is possible to optimize certain photoreactions by 
judicious choice of the irradiating energy. This technique 
is not possible with neutron activation analysis. 

Photon activation is the process involving nuclear reac- 
tions such as (y,n), (y,T), (&\np), et cetera. The 

high energy gamma rays are usually produced by the collision 
of a high energy electron beam with a high Z material such as 
tungsten, which results in a broad frequency distribution of 
photons, called the bremsstrahlung spectrum. Reported values 
in the literature for electron beam energy range from 6-115 
Mev and for beam current 5 uA to 1mA (See Table I). 

Table I. Survey of Electron Beam Energies and Currents 
Used in Photon Activation Analysis 



Principal 


Electron Beam Electron Beam 




Author 


Energy Mev 


Current uA 


Reference 


Kaminishi 


6 


1000 


21 


Chattopadhyay 


15-44 


25-225 


12 


Mulvey 


22 


250 


22 


Wilkens 


22 


30 or 45 


23 


Aras 


30 


50 


15 


Lutz 


35 


20 


20 


Hi slop 


35-40 


5 


1 


Ricci 


105-115 


210 


24 


One of the 


attractive 


features of photon 


activation 


analysis (PAA) 


is its use as 


a nondestructive technique for 


multi -element 


analysis. 


This aspect along 


with other 
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characteristics of the method are discussed by Hislop in his 
review: PAA of biological and environmental samples [Ref. 1: 
p. 1159]. Included in the review is a table listing elements 
found in the following materials: tobacco, tree bark, hair, 
kale, blood, soil, bone, air particulates, and urine. Brazil 
nuts [Ref. 2: p. 99], rubber [Ref. 3] and moonshine whiskey 
[Ref. 4] have also been investigated using photons. 

The focus of this study will be petroleum. The analysis 
of petroleum is motivated by a desire to use photon 
activation analysis as a technique to identify the source of 
an oil spill. It is known, for example, that the trace 
elements present in oil depend strongly on the source of the 
oil. The approach to be taken will be to first analyze 
environmental matrices such as charcoal and soil. The 
techniques developed are then extended to an exploratory 
analysis of petroleum. 
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II. THEORETICAL CONSIDERATIONS 



A. GENERAL TECHNIQUE 

The basic concept of either neutron or photon activation 
analysis is the production of radioactive nuclides as a 
result of some induced nuclear reaction. As the radioactive 
nuclide decays, selected radiations such as gamma rays are 
measured which will aid in the identification of the original 
elements in the sample. By judicious choice of the 
irradiation and counting variables it is possible to optimize 
the study of individual elements. 

There are five general steps in activation analysis: (1) 
the selection of an appropriate nuclear reaction; (2) the 
preparation of samples for irradiation; (3) the irradiation; 
(4) the post-irradiation assays; (5) the evaluation of the 
experimental data. The general principles for each of these 
steps, as applied to photon activation analysis, will be 
reviewed here. More detailed discussions will follow in 
subsequent sections. 

1. Nuclear reactions to be expected in photon activation 
are a function of the physical, chemical, and nuclear 
properties of the sample matrix and activation products. The 
two general situations are the investigation of a sample for 
known elements and the investigation of a sample for unknown 
elements. In either case the importance of interfering 
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and/or competing reactions must be considered. In selecting 
an optimum reaction a necessary condition is the availability 
of adequate sources of nuclear data. 

2. One of the advantages of photon activation is that it 
usually requires little or no sample preparation. 
Specifically, pre-irradiation chemical processing is not done 
except in special applications. Samples to be irradiated may 
be pure elements or a mixture of elements either in solid or 
liquid form. The usual containers for irradiation are 
standard polyethylene vials and aluminum capsules. 

3. The irradiation itself is accomplished with high 
energy photons. While there are a number of ways to produce 
photons the most common method is to use an electron linear 
accelerator. The variables of importance in the irradiation 
are the energy of the electron beam, the time of irradiation, 
the intensity of the irradiation, the number of target atoms 
and the geometry of the target. This method will be 
explained in more detail in Chapter II-C. 

In general a higher irradiating energy will give higher 
yields. However, competing reactions will also be enhanced. 
Therefore photon activation experiments are most often done 
with a lower irradiating energy (20-40 Mev electron beam) and 
higher beam intensities (20-250 uA beam current), as compared 
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to the typical 100 Mev, luA beam produced by the Naval 
Postgraduate School Linac. 

4. Post-irradiation assays are of two types: absolute 

or comparative. In the absolute assay accurate values of all 
variables (i.e., flux, cross sections, detector efficiency, 
et cetera) must be known. The comparative technique requires 
relative measurements only. In this method a sample and a 
standard are irradiated simultaneously. By taking the ratio 
of activities the need for certain nuclear data is 
eliminated. This gives more accurate results with more 
convenience than the absolute method. For this reason the 
comparative assay is used most often. 

With either method suitable counting equipment is 
required for the type of radiation to be observed. 
Frequently Nal(Tl) or Ge(Li) detectors are used to record the 
gamma ray spectrum of the decaying radionuclides. 

5. The final step in photon activation is the evaluation 
of the experimental data. For gamma ray spectroscopy this 
involves identifying the energy of the photopeaks observed in 
the gamma ray spectrum. When possible the half-life 
associated with a given photopeak should be measured. With 
this information photopeaks can be matched to particular 
radioisotopes and the original target nucleus deduced. 
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This final step should also include a discussion of the 
possible errors contributing to the accuracy and precision of 
the experiment. 

B. NUCLEAR REACTIONS 

A nuclear reaction occurs when a particle of sufficient 
energy interacts with a nucleus to induce a change in the 
nucleus. This process is different from the radioactive 
decay of a nucleus. A nuclear reaction is induced; 
radioactive decay is spontaneous. 

A nuclear reaction is written as 

A + a > B + b + Q 

where A is the target nucleus, a is the irradiating particle, 
B is the product nuclide, b is the particle emitted from A, 
and Q is the change in energy of the system. A shorthand 
notation for the above reaction is 

A (a , b) B 

The irradiating particles can be neutrons, electrons, 
photons or ions. For photon activation there are two 
parameters of fundamental importance. These are the energy 
of the incoming photon and the cross section for the 
reaction . 

The minimum photon energy for a reaction to occur is the 
threshhold energy. This is also called the separation energy 
by some authors since it is the energy at which a product 
nuclide begins to form. For (2f,n) reactions the threshholds 
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usually range from 7 to 18 Mev and they decrease with 
increasing atomic number [Ref. 5: p. 721], 

The reaction cross section is a measure of the probabi- 
lity a reaction will occur. It has units of barns where 1 
barn = 10 cm . Typical cross sections for photon reac- 
tions are on the order of millibarns. Like the threshhold 
energy, the cross section varies with the atomic number 
ranging from 1.6 millibarns for beryllium 9 to 1.58 barns for 
plutonium 239. Other cross section values can be found in 
the activation analysis handbook by R. C. Koch [Ref. 6]. 

The reaction cross section is also a function of the 
photon energy. For example, as the photon energy is increased 
in a (tf,n) reaction there is a greater chance that competing 
reactions will occur (i.e. ( 2T, p) , (tf,np), (}f,2n), et cetera. 
The result is a marked decrease in the cross section for 
(y,n). This phenomenon is known as the giant resonance and 
is pictured qualitatively in Figure 1. E T is the threshhold 
energy, the maximum cross section, and E 0 the resonance 
energy . 

Both the maximum cross section and the resonance energy 
are found to vary smoothly with the mass number from sodium 
to uranium. They are roughly proportional to NZ /A and 
A- 0 * 33 respectively [Ref. 7: p. 9]. A more exact expression 

is given by DeSoete [Ref. 5: 



for the resonance 
p. 728]. 



energy 
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E q = 40.7 A" 0,20 

The maximum cross section for 2 < Z < 20 varies between 1 and 
30 millibarns. For 20 < Z < 83 the peak cross section can be 
estimated by: 

(T (E 0 ) = 8.15Z - 156.4 mb 

Figure 2 displays the threshhold energy and resonance energy 
as a function of atomic number. It can be seen that the most 
efficient range of photon energies for activation purposes is 
15-25 Mev. 

For low Z nuclides the ( n ) and (8*,p) nuclides are of 
the same order of magnitude. However as Z increases the 
coulomb barrier of the nucleus hinders the emission of a 
proton {Ref. 5: p. 723]. Therefore the ( if , n) reaction is 
more probable for Z > 30. 

The nuclear reactions applicable to photon activation 
analysis are of two types. The first type are 
reactions which lead to isomeric states of an element. The 
threshhold for these reactions are usually well below 10 Mev. 
Because the excited nucleus decays back to the ground state 
by emitting a characteristic gamma, the reaction is 

attractive for both its selectivity and its non-destructive 
feature . 

The second type of reaction includes (2f,n), (y,2n), et 

cetera reactions where neutrons are ejected from the nucleus 
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cross section (mb) 



B 8 10 IS 14 16 18 80 22 24 26 28 30 32 34 36 38 




Gamma ray energy (Mev) 

Figure 1. Giant Dipole Resonance of Pb208 [Ref. 25: p. 736] 




o Threshhold energy 
x Maximum cross section energy 



Figure 2. Threshhold Energy and Energy of Maximum Cross 
Section for (Jf,n) Reactions [Ref. 8: p. 95] 
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and (2T,p), (Jf,T), (2f,np), et cetera, reactions where charged 
particles are emitted. These reactions are considered pri- 
mary reactions, that is, a nuclear reaction induced in an 
original sample element by the principle irradiating 
particle . 

Besides primary reactions, secondary and second order 
reactions are also possible in photon activation. A 
secondary reaction is a nuclear reaction in an original 
sample constituent by particles created during primary 
reactions. A common reaction would be a (n,^r) event. In this 
case neutrons produced from (}f,n) reactions may react with 
nuclei from the original sample. 

Second order reactions are of two types. The first type 
is one in which the principle irradiating particles interact 
with primary reaction products. Fortunately, this type of 
interference is not expected with photons since the reaction 
cross sections are very small [Ref. 6: p. 10]. The second 
type of second order reaction occurs when primary products 
undergo spontaneous decay. This type of interference may be 
important if neighboring elements of a matrix are being 
analyzed . 

To summarize nuclear reactions, it can be expected that 
with increasing photon energy a number of gamma reactions are 
possible. Which one will predominate is a function of photon 
energy, photon flux, reaction cross sections, abundance of 
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target nuclides and the length of irradiation. An excellent 
review of many reactions possible with photon activation 
analysis is given by Toms [Ref. 9, 10]. Figure 3 is a dis- 
play of possible primary gamma reactions. 

C. PRODUCTION OF PHOTONS 

There are four sources of photons for activation 
analysis: positron annihilation, nuclear reactions, 

radioactive isotopes and electron accelerators. Generally 
the least useful of these methods is positron annihilation, 
since the photon flux is about 10” 4 to 10” 8 times less than 
the bremsstrahlung flux available from electron accelerators 
[Ref. 11: p. 10]. 

The early work in photon activation analysis was 
accomplished with either natural or artificial radioisotopes. 
Borshov and co-workers, for example, in the late 1930's 
determined beryllium by prompt neutron emission with a radium 
source [Ref. 8: p. 93]. Other radioisotopes used to generate 
photons are Na24, Co60, and Sbl24. Strengths of the sources 
range from 0.2 to 30,000 curies. 

Production of photons by nuclear reactions is also 
possible. Del Bianko, for instance used monochromatic gamma 
rays from the H 8 (p,2f) He 4 reaction to study photonuclear 
cross sections [Ref. 7]. Two other possible reactions are 
Li 7 (p, <f) Be 8 and B 11 (p,^)C 12 . The gamma ray fluxes obtained 
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Figure 3. Possible Gamma Reactions 



from these reactions, however, are very small and the gamma 
energy is fixed [Ref, 5: p. 732]. 

The fourth method of producing photon fluxes is by 
colliding a high energy electron beam with a high atomic 
number material. The inelastic scattering of the energetic 
electrons by the coulomb attraction of the positively charged 
target nucleus leads to the emission of a continuous spectrum 
of photons. This process and the radiation produced is known 
as bremsstrahlung (braking rays). The maximum energy for the 
gamma rays is equal to the highest kinetic energy of the 
irradiating electrons. Besides radiative losses 
(bremsstrahlung), electrons passing through matter also lose 
energy by excitation and ionization. This process has two 
effects. One, it results in the heating of the target and/or 
sample. Two, the number of electrons converted to photons is 
reduced. For thick targets (length > 1.5 radiation lengths) 
the fraction of the electron energy converted to 
bremsstrahlung can be calculated from an expression given by 
Koch and Motz : 

£ = 1 - In (1 + 1.2 x 10" 3 EZ)/1.2 x 10" 3 EZ 
where £ is the fractional efficiency, A is the atomic number 
and E is the electron energy in Mev [Ref. 8: p. 94]. For 
lead (Z = 82) and E = 100 Mev electron beam, £ = 75%. 

An important feature of bremsstrahlung is its anisotropic 
character as the produced gamma rays are peaked forward 
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[Ref. 5: p. 732]. The number of photons in 

direction increases with increasing electron 
angular width is roughly l/tf in radians, 

( (electron) /me 4 ) . 



the forward 
energy. The 
where Y ~ E 




Figure 4. Schematic Sketch to Illustrate the Spatial 
Distribution of Bremsstrahlung 



Van de graaff generators, synchrotrons, betatrons, 
microtrons, and linear accelerators are all capable of 
producing electron beams for use in bremsstrahlung. The most 
commonly used apparatus for photon activation, however, is 
the linear accelerator. This is because a requirement for 
trace analysis is a high intensity beam. Cha t topadhyay and 
Jervis recommend a photon flux of about 10^ 3 to lO^- 4 
photons/cm 2 /sec. [Ref. 12: p. 1639]. DeSoete, et. al. 
suggest that an integrated current of more than 10 uA and 
energies up to 40 Mev be used for activation analysis 
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[Ref. 5: p. 278]. These levels of intensities and energies 
are easily obtained with most linear accelerators. 

D. ACTIVATION EQUATION 

The activation equation in its simplest form can be 
expressed as 

(1) A« w 

which states that the activity of an element is proportional 
to its weight. In more detail equation (1) is written as 

(2) A (o) = NC“4>(l-e" XT ) 

where A(o) = the activity of a radionuclide at the end of an 
irradiation 

N = the number of target atoms 
<r* = the reaction cross section 

4> = the flux of irradiating particles 

A = the decay constant 

T = the irradiation time 

For a particular radionuclide equation (2) is modified to 
account for the fractional abundance, f, of the isotope being 
analyzed . 

(3) A (o) = _Wf<T<j> (l-e _AT ) N q 
M 

Here the number of target atoms has been expressed as W/M N Q 
where W is the weight of the element present in the sample, M 
is its atomic weight, and N Q is Avogadro's number. 

Usually the activity will not be determined until some 
time t after the irradiation. Therefore another factor must 
be added to equation (3) to account for the radioactive decay 
of the isotope being measured. The activity is then 
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(4) A ( t) =_Wf <rcb(l-e~ AT ) e” At N Q 

M 

This equation is rearranged to yield 

(5) w = MA (tie At 

f<r4>(l-e" AT )N 0 

In theory if all the variables on the right hand side of 
equation (5) are known then the weight of an element can be 
calculated directly. This is an example of an absolute 
assay. In practice there is enough uncertainty in the 
reaction cross section and flux that in most cases a 
comparative technique is used. 

As explained previously, the comparative assay involves 
the simultaneous irradiation of a sample and a standard. 
Using equation (5) the weights of a desired element in a 
sample X and a standard S are 

< 6a > w x * - 

frp(l-e )N Q 

( 6 b) « s - MA 

f<r9(l-e Ai )N Q 

Assuming a homogeneous irradiation of the sample and 

standard the ratio of W x /W s gives 

t: 

Hrs- 



(7) W v - A y (t x)e l~ x 



W 3 A s ( ts ) e 

The activities in equation (7) are determined 
experimentally by counting the sample and standard under 
comparable conditions. W s is known, so W x can then be 
calculated. 
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1 . Sensitivity Calculation 



A modified equation (5) can be used to measure the 
sensitivity of an element for activation analysis. The 



minimum weight of an element detectable in an irradiated 



sample is 



(8) W, 



mm 



= MAm f t) e 

-f<r4>(l-e~' AT )N 0 YF£ 



where M = atomic weight of the element 

Am(t) = minimum detectable counting weight at the time 
of counting 

e A t = decay from end of irradiation to start of 
counting 

f = isotopic abundance 

<T = reaction cross section 

$ = irradiation flux 

1-e - '*™ = saturation factor for the length of the 
irradiation 

Y = chemical yield if radiation chemistry is used 

F = self absorption losses 

£ = overall counting efficiency 



A few of these variables merit additional comments. Am(t) is 



the minimum resolvable counting rate above background and is 
often assumed to be 2B where B is the background [Ref. 13: 
p. 242]. 

The overall counting efficiency is given by 

£ = ci 

D 



where Ci is the number of counts of radiation i and D is the 
number of disintegrations of the nuclide in the same time 
interval [Ref. 13: p. 243 ]. It depends both on the 

radiation-energy efficiency of the detector and on the decay 
scheme of the radioisotope being measured. For example a 
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Nal(Tl) detector may have a 20% efficiency for 2.13 Mev gamma 
rays, but the radioisotope Cl34m decays only 40% of the time 
with the emission of 2.13 Mev gamma rays. The overall 
counting efficiency is then 

£ = 0.20 x 0.40 = 0.08 

The counting geometry will affect the overall counting 
efficiency. Usually it is determined from a calibration 
curve of £ vs gamma ray energy for a fixed set of counting 
conditions . 

From equation (8) it can be seen that the maximum 
sensitivity (i.e., minimum weight) is obtained when the 
variables in the numerator are minima and the variables in 
the denominator are maxima. As an example of the 
sensitivity calculation, an activation analysis for the 
determination of lead in petroleum will be evaluated. The 

0.899 Mev photopeak from the Pb206 Of, 2n) Pb204m reaction is 
used for analysis. 

1. Counting system: 

Nal(Tl) detector, total B (photopeak) = 17 cpm 
Am(t) = 34 cpm = .568 cps 
overall counting efficiency = 32% 

2. Post-irradiation processing: 
elapsed time to counting : 3 hours 
Y = 1 (no chemical processing) . 

F = 1 (no self absorption) . 

3. Irradiation conditions: 

<r = 85 mb 

$ = 2 x 10 J ' J photons/cm ^sec . 
irradiation time = 1 hour 
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4. Sample data: 

isotopic abundance of Pb206 = .236 
sample size = 1 gram oil 

assume no major interfering constituents, 
atomic weight of lead = 207.2 g/mole. 

/\= .0108 min" 1 

Equation (8) yields: 



w min = 3Q7. 2 f . 568) . fi- . 0103(1 fQ J nin^/ ' g TH T* 

.236 (85x10 (2xlO iJ ) (1-e" * 0103 {60) ) (6.02xl0 23 ) (.32) 



W, 



mm 



= 2.1 x 10" 3 grams 



Usually it is the concentration of an element in a sample 
which is desired rather than the absolute amount in the 
sample. This concentration is obtained by dividing W m ^ n by 
the sample weight and multiplying by 10° for ppm and 10^ for 
percent. 



< 9a > c min = *min x 10 6 ppm 
S 

< 9b > ^in - «min * 11)2 * 

s 



For the above calculation the concentration of lead which 
could be determined in a sample under the stated conditions 
would be .021 ppm. Although this is an idealized calculation 
because no interfering reactions were assumed, the excellent 
sensitivity calculated is not unreasonable for activation 
analysis. Sensitivities of nanograras are frequently reported 
in the literature for activation analysis. 



E. GAMMA RAY SPECTROSCOPY 

The most common method for collecting data in photon 
activation analysis is to use either a Nal(Tl) or Ge(Li) 
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detector coupled to a multichannel analyzer. To evaluate the 
gamma ray spectrum obtained it is helpful to understand the 
basic principles of gamma ray spectroscopy. This section 
will be a review of those principles. 

Not considering cost/ the choice of using either a 
Nal(Tl) or Ge(Li) detector is driven by the type of analysis 
to be done. Both detectors have advantages. For example, 
the Ge(Li) crystal has an intrinsic high resolution. The 
spectrum obtained is characterized by very narrow peaks 
(typically less than 1/10 the width of Nal peaks) that lie on 
relatively simple continuous compton distributions [Ref. 14: 
p. 636]. The peaks are thus widely separated and completely 
non-interfering. A Na24 spectrum obtained with a Nal(Tl) and 

a Ge(Li) detector is compared in Figure 5. However, the 

/ 

efficiency of the Ge(Li) crystal is only a small fraction of 
the Nal(Tl) crystal. A representative number would be 6-10% 
compared to a standard Nal(Tl) crystal. Thus the choice of a 
detector is a compromise between resolution and efficiency. 
For low level counting of a single element in a relatively 
simple matrix, high efficiency is needed and hence Nal(Tl) 
would be the preferred choice. For multielement analysis of 
complex environmental or biological samples resolution is the 
most important criterion and Ge(Li) detectors must be used. 
Regardless of which detector is used, the general operating 
principle is the same. 
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Figure 5. Comparison of Na24 Spectra Taken with Nal(Tl) and 
Ge(Li) Detectors [Ref. 26*] 
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Gamma rays interact in a crystal (Nal(Tl) or Ge(Li)) in 
essentially three ways: (1) the photoelectric effect/ (2) 

the compton effect/ and (3) pair production. All three of 
these interactions produce secondary electrons which deposit 
their energy in the crystal/ resulting in a scintillation in 
the Nal(Tl) crystal and an electronic pulse in the Ge(Li) 
crystal. With Nal (Tl) , the crystal is coupled via a 
photocathode to a photomultiplier tube to convert the 
scintillations into electronic pulses. With both crystals/ 
the electronic pulses are amplified, sorted by amplitude, and 
stored in a pulse-height analyzer. Since the pulses are 
proportional to the incident energy deposited in the crystal, 
the end result is a gamma ray spectrum. If more than a few 
radionuclides are decaying in a sample, the resulting 
spectrum can be very complex. 

The major characteristics of a gamma ray spectrum are: 

a. photopeaks 

b. the compton edge and distribution 

c. backscatter peaks 

d. sum peaks 

e. single and double escape peaks 

f. tube noise. 

These peaks are illustrated in Figure 6. Their origins can 
best be explained by discussing Figure 7. 

Gamma rays with energies up to about 0.1 Mev interact 
within a Nal(Tl) crystal predominantly by the photoelectric 
effect. The photoelectric effect is a direct interaction of 
photons with the orbital electrons in an atom. Up to about 
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Figure 6. Hypothetical Gamma Riy Spectrum 




Figure 7. Interactions of Gamma Rays with a Nal(Tl) Crystal 
[Ref. 14: p. 640j 
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so in Nal (Tl) , the 



0.5 Mev this effect is proportional to Z^, 
photons interact mostly with the iodine atoms. This is an 
inelastic collision the result being the ejection of an 
electron from the atom with kinetic energy. 

(1) E e = Ey - B e 

where B e is the binding energy of the ejected electron in its 
orbital shell. The iodine atom will be left with a vacancy 
in its K shell which when filled by cascading electrons gives 
rise to a characteristic x-ray of 0.028 Mev. 

The energy of the photoelectron in equation (1) is 
quickly deposited in the Nal(Tl) crystal. The K x-ray may or 
may not deposit its energy in the crystal depending on where 
the interaction took place. Both possibilities are depicted 
by ^ and %2 i- n Figure 7. If the interaction occurs near the 
surface the K x-ray may escape the crystal. If the interac- 
tion occurs elsewhere in the crystal the K x-ray is almost 
always captured. Thus, at low gamma ray energies two peaks 
can occur in the gamma ray spectrum: one at the incident 
gamma ray energy and one at 0.028 Mev below this. Because 
the difference between these two peaks is small in absolute 
magnitude, these two peaks are not observed separately at 
gamma ray energies above 0.170 Mev [Ref. 14: p. 641]. Only 
one peak would then occur at a pulse height corresponding to 
the incident gamma ray energy. This total absorption of the 
incident gamma ray represents the photopeak. 
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The next important process contributing to the formation 
of a gamma ray spectrum is the compton effect. This effect 
is similar to the photoelectric effect. However/ an elastic 
rather than an inelastic collision occurs between the 
incident gamma ray and an electron. A compton electron will 
be the result from this collision and the gamma ray will 
emerge from the collision in a new direction and with reduced 
energy. This degraded gamma ray can then undergo additional 
compton scatterings or interact by the photoelectric effect. 
This process is illustrated by ^3 in Figure 7. The total 
energy deposited in the crystal due to ^ would be e c + e c ' + 
ep + x. The total time for these interactions is so fast 
compared to the decay time of a light pulse from the Nal(Tl) 
detector that only one light pulse corresponding to the 
photopeak will be emitted by the crystal [Ref. 14: p. 642]. 

By solving the energy and momentum equations for the 
compton effect/ the kinetic energy of the compton electron/ 
T, and the energy of the scattered photon, E’, can be 
expressed as 

(2a) T = E q 

(i-cose) Eq 

(2b) E 1 = E„ 

1 + (l-Cos^ E~ 

mo ^ 2 

9 

where E Q is the energy of the incoming photon, m 0 c is the 
rest energy of the electron, and -O-is the angle through which 
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the photon E' is scattered relative to the original direction 

of the photon E Q . From equation (2a) it is evident that as 

E^ approaches 0, T approaches 0. That is to say, there is 

no minimum kinetic energy for the compton electron. However 

a maximum exists at ^ = 180°. 

(3a) Tniax = E n 

TT> m n c z 

The corresponding photon energy for 180° is given by 

(3b) _E^_ 



m o c 

Referring back to the case of ^3 in Figure 7, it is 
possible for the scattered photon 83 ' or # 3 " to escape from 
the crystal. If this situation occurs an energy of E Q - E' 
is deposited in the crystal. Since E' can vary from E 0 when 
-©*= 0 to E Q - T max ( w ^en *©■= 180°) a continuous distribution 
of peaks will appear in the gamma ray spectrum from 0 to 
T max- T max i s ca Hed the compton edge. 

The final mechanism important in gamma interactions is 
pair production. This event usually occurs in the coulomb 
field of the nucleus and in the process the gamma ray 
disappears and an electron positron pair is created. The 
minimum gamma ray energy for pair production to occur is 1.02 
Mev (twice the rest mass energy of an electron). Pair 
production processes become significant in Nal(Tl) when the 
initial gamma ray energy is greater than 2 Mev. 
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V 4 in Figure 7 illustrates the possible events which can 
occur in pair production. The energy of the electron in the 
electron-positron pair is quickly deposited in the crystal. 
The positron readily annihilates with an electron. This 
interaction produces two photons of .511 Mev traveling at an 
angle of 180° to one another in an arbitrary reference frame. 
These two photons may then interact by the compton or 
photoelectric process, or one or both may escape from the 
crystal. If neither photon escapes then a photopeak will 
appear in the gamma ray spectrum corresponding to the energy 
of If one .511 Mev photon escapes a peak will form at 
the energy ^4 - .511 Mev. This peak is called the single 
escape peak. If both .511 Mev photons escape from the 
crystal then a double escape peak forms at -1.02 Mev. 
Continuous compton distributions are also present for each 
peak . 

Another peak possible in a gamma spectrum is the sum 
peak. This occurs when the detector counts two gamma rays as 
one. These peaks become noticeable if the counting rate is 
high. 

Besides the photopeaks, compton distributions, single and 
double escape peaks, and sum peaks produced by direct 
interactions of gamma rays from a source, secondary effects 
due to the inevitable surroundings tend to distort a gamma 
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ray spectrum. The various events possible are shown in 
Figure 8. 

Compton scattering of the primary radiation from the 
surrounding material contributes to the pulse height 
distribution from zero amplitude to a maximum depending on 
the angle of scattering. The most common effect is 
backscatter ing (i.e. 180° scattering) behind the source. The 
theoretical location of the backscatter peak is 

e BS e o T max 

where E BS is the energy of the backscatter peak and E Q and 
T max are the incident gamma ray energy and compton edge 
respectively. 

Annihilation radiation at .511 Mev is also common in a 
spectrum especially if a source decays by positron emissions. 
Pair production in the shielding material will also produce 
annihilation radiation. 

The last probable interference from surrounding materials 
are x-rays produced by the photoelectric effect. A .074 Mev 
peak corresponding to the Pb x-ray is common. 

Finally, tube noise from thermionic emissions and natural 
background from cosmic rays and surrounding radioactive 
materials complete the contributions to a gamma ray spectrum. 

In summary, a spectrum typically consists of the 
following: photopeaks, compton distributions, single and 

double escape peaks, sum peaks, backscatter peaks, and noise 
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Figure 8. Summary of Nuclear Events in the Vicinity of a Nal(Tl) Detector (Ref. 27: 



effects. 



Table II and Figure 6 summarize the above 



discussion. 



Table 


II. Origin of Peaks 


in a Gamma 


Ray Spectrum 


Photon 

Escape 

Energy 

(Mev) 


Origin 


Spectrum 

Peak 

Energy 

(Mev) 


Spectrum 

Peak 

Name 


0 


Total absorption 


E o 


Photopeak 


0.511 


Pair production 
and escape of one 
.511 Mev photon 


E q -0 . 511 


Single escape 
peak 


1.02 


Pair production 
and escape of two 
.511 Mev photons 


E q -1 . 02 


Double escape 
peak 


E min 


Compton 180° 
scattering 


T max 


Compton edge 


El min to E 


0 Compton 
scattering 


0 to T max 


Compton 

distribution 


E Q with 
E BS r ?" 


External Compton 
180° scattering 


e bs 


Backscatter 

peak 



entering 



Adapted from J. B. 3 i r k s , T he T he ory and Pr a ct i ce of 
Scintillation Counting , Macmillan, 1964, p. 470. 

F. PRECISION AND ACCURACY 

The emphasis on precision and accuracy in the analysis of 
experimental data depends on the type of analysis to be done. 
For example, a very accurate result was required by Guinn and 
co-workers for their determination of lead in moonshine 
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whiskey for forensic purposes [Ref. 4]. However, in their 
exploratory studies of biological samples, drugs, glass, 
soil, and paper, the main effort was required for the 
identification of the elements found in the samples. Thus a 
rigorous error analysis was not required, implying less 
overall accuracy for the procedure. 

In general, the need for precision and accuracy can be 
divided into four categories, according to the question to be 
answered [Ref. 13: p. 240], 

1. Qualitative: is element 2 present in the sample? 

2. Threshhold: is element Z present in amounts greater 
than some given amount? 

3. Relative: is element Z present in small or large 
amounts? 

4. Quantitative: exactly how much of element Z is 
present in the sample? The greatest need for precision and 
accuracy is when quantitative results are required. 

There are a number of errors in photon activation 
analysis which affect the overall accuracy. The main two 
are: (1) changes in the irradiation conditions, and (2) 
interfering nuclear reactions. These errors and other 
sources of errors are summarized in Table III. 

Changes in the irradiation conditions refers to flux 
gradients in the sample. Careful attention to the sample 
geometry and close monitoring of the irradiating beam energy 
will minimize this type of error. 
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Table III 



. Sources of Error in Photon Activation Analysis^- 







Source 


Estimated Error 


1. 


Chemical 


% 




a . 


Sample weight 


-kl 




b. 


Standard weight 


±2 




c . 


Yield (if seperation used) 


±2 


2. 


Irradiation 






a . 


Irradiation time (< 1 min) 


±3 




b. 


Nonhomogeneity of photon flux 


±2 




c . 


Electron beam energy 


±1 


3. 


Nuclear data 






a . 


Half-life 


±2-10 




b. 


Decay scheme 


±2-50 




c . 


Cross section 


±5-30 


4. 


Nuclear reactions 






a • 


Competing 


Variable 




b. 


Interfering 


Variable 


5. 


Counting 






a . 


Detector calibration 


±3 




b. 


Counting rate < 10 3 dps 


±4 




c • 


Geometrical factors 


±1 



1 Adapted from J. P. Call, et. al., "The Accuracy of 
Radioactivation Analysis", in Modern Trends in Activation 
Analysis (Texas A & M University, College Station, 1965), 
pp. 258-258 . . 



Interfering nuclear reactions occur when the same 
radionuclide is produced by more than one mechanism. An 
example is the production of Na24 by the following reactions. 

Mg 2 5 ( Y^p) Na24 
A127 ( )f,He 3 ) Na24 
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This type of interference can be reduced or eliminated by the 
proper choice of irradiating energy. For example, the 
threshhold energy for the (}f,He^) reaction for A127 is 23.7 
Mev. Thus by using 20 Mev electrons to produce a 
bremsstrahlung flux of photons, no Na24 will be produced from 
A127 . 

Another type of interference probable in photon 
activation analysis is overlapping gamma rays in the gamma 
ray spectrum. Various techniques are used when this type of 
interference occurs. If, for example, the interference is 
due to a short lived nuclide, the contribution of the short 
lived nuclide to the photopeak is minimized (if not 
eliminated) by waiting five half-lives before counting. On 
the other hand, if a long lived nuclide interferes, it can be 
produced in smaller quantity by limiting the length of the 
irradiation. 

When all possible sources of errors are minimized, the 
overall accuracy of photon activation analysis is around 
± 10%. A value of ±10 to ±15% is given by Aras, et. al. for 
their analysis of atmospheric particulate material [Ref. 15: 
p. 1407]. Chattopadhyay and Jervis estimated that the total 
error involved in analysis can be restricted to ±7 to ±8% in 
some instances, and to ±10% in all cases [Ref. 12: p. 1638]. 
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III. EXPERIMENTAL PROCEDURES 



A. GENERAL 

The experiments done for this thesis can be classified in 
one of two groups. The first group involves experiments done 
with petroleum samples. The second group includes other 
environmental samples such as charcoal, charcoal residue, 
aluminum foil, soil, and a crystal of bismuth germanate 
(BGO) . 

Certain experimental conditions were constant for all 
experiments. In all cases except BGO, samples were 
irradiated with a bremsstrahlung flux of photons produced by 
the collision of a high energy electron beam with a 5.3mm 
thick lead target. The electron beam was generated by the 
Naval Postgraduate School's electron linear accelerator. The 
Linac is capable of producing electron energies between 20 
and 100 Mev. The typical beam intensity was 5 x 10^ 
electrons/second . 

Gamma rays emitted from radionuclides produced in the 
irradiation were counted with a 3 inch x 3 inch hermetically 
sealed Harshaw Nal(Tl) detector. The FWHM range for the 
0.667 Mev gamma of Cs-137 was measured as .060 Mev, implying 
9% resolution. For the 1.332 Mev gamma of Co-60 the FWHM 
value was .084 Mev for 6.3% resolution. 
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Pulses from the Nal(Tl) detector were amplified with an 
ORTEC 485 amplifier and then sorted and stored by a TRACOR 
NORTHERN 7200 multichannel analyzer. A two point energy 
calibration program of the multichannel analyzer made it easy 
to identify the gamma ray energies corresponding to all 
photopeaks. When possible the half-life associated with a 
given photopeak was measured. 

B. COLLECTION OF SAMPLES 

Samples were collected from a variety of sources. The 
charcoal sample was obtained by crushing a briquette of 
Kingsford "Match Light". The charcoal residue was also 

A 

Kingsford's although after the barbecue. Two types of 
aluminum foil were used: Safeway and Reynolds Wrap. The 
soil sample was a commercial potting soil purchased locally. 
The BGO crystal analyzed was a fragment of a larger BGO 
crystal obtained from Harshaw Chemical Corporation. Finally, 
the oil samples came from one of three sources. 
Approximately 50 grams of used engine oil from a fried's car 
was contributed for analysis. Another 34 samples consisting 
of 12 crude oils, 6 lube oils, 6 #6 fuel oils, 3 bunker fuel 
oils, 6 #2 fuel oils, and 1 #5 fuel oil were provided by the 
U.S. Coast Guard's Central Oil Identification Laboratory. 
Five crude oils were also available for analysis from the 
U.S. Coast Guard Research and Development Center, Groton, 
Connecticut. 
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c 



IPAA -GENERAL SAMPLES. 



1. Aluminum Foil 

As a check on the homogeneity of an irradiation, it 
is advisable to use a flux monitor. A suitable material to 
use in major and minor element analysis is household aluminum 
foil. Two brands, Safeway and Reynolds Wrap, were easily 
obtainable. An evaluation of the gamma ray spectrum obtained 
after similar irradiations indicated no significant 
difference between the two types of foil. Indeed, the only 
real difference was about twenty cents in price. 

The two possible photoreactions of aluminum are 
listed below. 



Principal 





Photo 

reaction 


Residual 


Gamma 


Rays 


Target 


Sip 


Nucleus 


Energy 




(abundance) 


Mev 


1/2 Life 


Mev 


Percent 


13 A127 


He 3 


Na24 


1.368 


99.0 


100.0 


23.71 


15. Oh 


2.754 


99.0 




2p 


Na25 


0.98 


15.0 




22.42 


59s 


0.58 


14.0 








0.40 


14.0 








1.61 


6.0 


To use the 


foil as a flux 


monitor , 


the sample 


vials are 



wrapped with a sheet of aluminum foil. At a suitable time 
(around 4-5 hours) after the irradiation, the number of 
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counts in the 1.36 Mev peak due to Na24 are recorded. The 
activity at the end of the irradiation is then calculated by 
the following expression. 



(1) A(o) 



A(t)e 



At 



dps 



where A(o) = activity at the end of the irradiation 

A ( t) = activity at time t 

A = decay constant 

£ = total detector efficiency 

The flux is then calculated by 

(2) <j) = A.L o) . photons 

N 0 <r(l-e -/Vr ) f cm^ .sec 

where N Q = number of target atoms 

<r* = reaction cross section 

- saturation factor 
T = time of irradiation 

f = isotopic abundance 

It is important to realize that the flux given by 
equation (2) is not the actual flux present in the sample. 
This is because both the reaction cross section and the 
photon flux are functions of the photon energy. Strictly 
speaking, the quantity which should be solved for in equation 
(2) is the product . To be correct, the production rate 
for a photonuclear reaction must be calculated by summing the 
product of the cross section and the photon flux in an energy 
increment over the energy range of interest [Ref. 16: 
p. 426]. However, as a quick check on the homogeneity of an 
irradiation, equation (2) is still useful. As an example, 
consider two separate irradiations of a sample. Using the 
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procedure described above, and <$> 2 are calculated for the 
two aluminum foils. The ratio of is given by 

(3) = At (o) Noq 

”(J>2 A 2 (o) No ^ 

Equation (3) is useful as an order of magnitude check 
for an analysis. In one of the experiments for oil, 4*i was 
calculated as 4.4 x 10 3 photons/cm 3 /sec. 2 was calculated 
as 3.3 x 10® photons/cm^/sec. The ratio was thus 1.33 
indicating a 33% difference in flux. This variance of flux 
is common in beam experiments and is another reason the 
comparative assay is favored over separate irradiations. 

Another use for the aluminum foil is as a check for 
thermal neutron interfering reactions. For example, the 
aluminum foil from a one-half hour irradiation of a sample 
with 100 Mev bremsstrahlung radiation was counted with the 
Nal(Tl) detector. The following photopeaks were observed. 

Peak Energy (Mev) Production Mechanism 

0. 511 Annihilation Peak 

0 . 

1 . 

2.74 } A127(^,He 3 ) Na24 

Mg27 is produced in the following manner. 

Primary reaction A127(y,p) Mg26. 

Subsequent reaction Mg26(n,^0 Mg27. 

The implication of the above is that for major and 
perhaps minor constituents, products other than those 



8481 

015} Mg26(n,JO Mg27 
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produced by photoreactions may be encountered in a spectrum. 
This, of course, complicates the analysis. 

2. Charcoal Residue 

A 25 milliliter sample of charcoal residue was 
collected from the remains of a barbecue. The sample was 
irradiated for 10 minutes using a 100 Mev electron beam and a 
5.3mm lead convertor. Gamma ray spectrums were recorded over 
a 45 day period at various times after the irradiation. The 
nuclear reactions observed are listed in Table IV. Some 
insight into the thought processes involved in activation 
analysis will be provided by a discussion of how the elements 
in the matrix were deduced from the experimental data. 

The first clue as to what elements are in a sample 
comes from either previous knowledge of the sample or a 
literature search. For charcoal residue it can be assumed 
that the matrix is similar to fly ash. Table V lists 
concentrations of various elements found in fly ash. It can 
be seen that the fly ash matrix is composed mainly of iron, 
silicon, calcium, potassium and magnesium. Thus, it can be 
expected that these elements will be found in the charcoal 
residue. Assuming these elements are in the charcoal residue 
the probable reactions are listed. (Toms [Ref. 9, 10] is 
useful for this task.) The things to consider in deciding 
which reactions are probable include the following: isotopic 
abundance, cross sections, half-lives of products. 
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TablelV. Reactions Observed in the Photoactivation of 
Charcoal Residue 



Target 

Nucleus 

(abundance) 


Photo- 

reaction 


Residual 
Nucleus 
1/2 Life 


Principal 
Gamma Rays 

Energy 

(Mev) Percent 


19K39 


( )Trn) 


K38 


2.170 


100.0 


93.08 




7.7m 






20Ca40 


( Y,np) 


K38 


2.170 


100.0 


96.97 




7.7m 






20Ca44 


( 4p) 


K43 


0.374 


84.7 


2.06 




22. 4h 


0.619 


82.4 








1.02 


2.0 


20Ca48 


( Vf n) 


Ca47 


1.297 


75.0 


0.18 




4 . 5d 






20Ca47 


beta minus 


Sc47 


0.159 


70.0 




decay to 


3 . 4d 






12Mg25 




Na24 


1.368 


99.0 


10.00 




15. Oh 


2.754 


99.0 


12Mg26 


U,np) 


Na24 


1.368 


99.0 


11.01 




15. Oh 


2.754 


99.0 


12Mg24 


OT,2n) 


Mg 22 


no gamma 


rays 


78.99 




4. OS 


observed 




Mg 2 2 


beta plus 


Na22 


0.511 


180.0 




decay to 


2. 6y 


1.275 


100.0 


26Fe56 


( ^np) 


Mn54 


0.834 


100.0 


91.66 




312. 5d 






26Fe57 


(V/P) 


Mn56 


0.847 


99.0 


2.19 




2. 6h 






14Si29 


U,p) 


A128 


1.78 


100.0 


4.70 




2.3m 






14Si30 


u,p) 


A129 


1.28 


93.0 


3.09 




6 . 6m 
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Table V. Trace Element Concentrations in Coal, Fly Ash, and Fuel 
Oil [Ref. 19: p. 244] 





Coal, 


Fly aah. 


Residual fuel 


oe akmwnt 


ppm 


Ppm 6 


oil, ppm 


Hg 


0.02-2 


0 . 1- <18 


0.002-0.4 


Be 


<0.1-0. 4 


1-7 


0.0005- <0 .5 


Cd 


0.7-<30 


2- <100 


0.003-1 


As 


0.25- <100 


2.8- <200 


0 . 2- < 1 


V 


5.5-10 


180-2000 


40-113 


Mn 


1.9-20 


150-500 


0.21-1 


Ni 


4- <40 


45-300 


20-90 


Sb 


0.04-00 


5.6- <100 


0 .003- <0.5 


Cr 


3.4-00 


80-500 


0.7-4 


Zn 


5— < 100 


70-1000 


0. 4-2.0 


Cu 


<0.4-10 


33-300 


0.2-1 


Pb 


1.8-00 


95-440 


1-4 


Se 


0.1- <15 


0.77-40 


0.02-0.15 


B 


5-15 


190-500 


0.002-0.2 


F 


<2-60 


<10-100 


0 004* 


Li 


0.3-000 


20-300 


0.02- <3 


Ag 


<0.1- <2 


0 .04- <3 


0.0006-0.1 


Sn 


0.19- <30 


1.9- <100 


0.01-5 


Fe 


1800-8000 


5.3%-26% 


10-20 


Sr 


46-160 


69- <1000 


<0 4-<0.5 


Na 


100-870 


500-6600 


<0.4-30 


K 


20-2200 


0 .5%-3 .1% 


0.8-5 


Ca 


5500-1000 


1.3%-5% 


7- <400 


Si 


3000-20000 


>10% -20% 


8-30 


Mg 


600-2000 


2200-44000 


2-3 


Ba 


<2-500 


110-700 


0.3-5 


P 


NA- 


NA 


NA 


S 


NA 


NA 


NA 
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irradiation energy, 



irradiation time, irradiation energy, and counting 
conditions. In this sense, probably refers to the 
probability that a gamma ray corresponding to a particular 
radionuclide will be observed. For example, the (Y,n) 
reaction of Mg24 produces 12.1 second Mg23 with a gamma ray 
at 0.439 Mev. while the probability of this reaction 
occurring is high, the probability of observing the 0.439 Mev 
gamma ray is small unless the sample is counted immediately 
after the irradiation. 

Once the probable reactions are listed, the sample is 
irradiated. Gamma ray spectra are then recorded at various 
times following the irradiation. Gamma ray energies are 
assigned to each photopeak and if possible half-lives 
corresponding to a given photopeak are measured. Using the 
gamma ray energy, the half-life, if measured, and the 
predicted reactions, the elements in a matrix can be deduced. 
In the case of the charcoal residue, the 0.347, 0.619 and 

1.02 Mev gammas positively identified K43, from which the 
presence of Ca44 could be deduced. Likewise, the presence of 
two prominent peaks at 1.37 and 2.75 Mev identified Na24. 
The most likely reaction producing Na24 is the (JT,p) reaction 
of Mg25. 



Following the above procedure, the majority of the 
elements in a sample can be easily determined. However, in 
the case of unexpected gamma rays, a bit of nuclear sleuthing 
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may be necessary. For instance, in specta taken of the coal 
residue sample three weeks after the irradiation, the 
following peaks appeared consistently: 0.511, 0.847, 1.27 
and 1.81 Mev. Careful analysis seemed to indicate that the 
peaks were not decaying, so no half-life information was 
obtainable except to say that the half-life must be greater 
than 30 days. Checking a list of possible photoreactions 
which would produce a gamma ray energy corresponding to any 
of the observed gamma rays yielded the (ft^np) reaction of 
Fe26 producing Mn54 with a 312.5 day half-life and a 0.834 
gamma ray. The observed 0.847 Mev gamma ray does not quite 
agree with 0.834 Mev but given the resolution of the Nal(Tl) 
detector used, a difference of 0.013 Mev is reasonable. 
However, in a 10 minute irradiation it does not seem 
plausible that much activity can be expected from a nuclide 
with a 312.5 day half-life. Nonetheless, it is concluded 
that the reaction did occur for four reasons. One, the 
presence of iron in the sample was identified by other 
reactions. Two, the intensity times abundance for the 0.834 
gamma ray is quite high (91.66%). Three, the reaction was 
observed by Kenneth Murray in a 50 mg sample of Hawaiian 
pumice with a 27 Mev bremsstrahlung irradiation for 5 hours 
[Ref. 10: p. 1]. Four, the reaction was observed in an oil 
sample spiked with 5000 ppm of iron. (see page 83). 
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The identification of the 0.511/ 1.27 and 1.81 Mev 
gamma rays proved more difficult. First, it was not possible 
to find a photoreaction which would produce a nuclide with 
the required gamma rays. A search was then conducted for 
possible (n,lf) reactions which might have accounted for the 
observed activity. This effort was futile. Finally, by 
accident, it was noticed that the gamma ray combination was 
identical to the decay scheme of Na22. (Note, 1.81 is a sum 
peak.) Na22 can be produced via the Mg24 (X',n)Mg23 reaction. 
Mg23 is not observed because its half-life is 4.0 seconds. 
However, the decay of Mg23 to Na22 will yield the sought 
after gamma energies. 

The above discussion again points out the degree of 
complexity which is introduced in a gamma ray spectrum if 
more than a few nuclides are decaying in a sample. It should 
be remembered that a proper choice of irradiation time, 
irradiation energy, and counting schedule can help simplify 
the analysis. 

3 . Charcoal 

The charcoal sample was obtained by crushing a 
briquette of Kingsford "Match Light". Enough charcoal was 
crushed to fill a 25 milliliter polyethylene vial identical 
to the type used for the charcoal residue sample. The sample 
was irradiated for 10 minutes using a 70 Mev electron beam 
and a 5.3mm lead convertor. 
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Typical elements known to exist in coal are listed in 
Table V. The major constituent is carbon. Minor elements 
include iron, potassium, calcium, silicon and magnesium. 
Thus, the charcoal matrix is very similar to the charcoal 
residue matrix. Comparing Table IV to Table VI, it can be 
seen that essentially the same reactions occur in both 
samples. 



Table VI. Reactions Observed in the Photoactivation 
of Charcoal 



Target 

Nucleus 

(abundance) 


Photo- 

reaction 


Residual 
Nucleus 
1/2 Life 


Principal 
Gamma Rays 

Energy 

(Mev) Percent 


6C12 


(J'/n ) 


Be 7 


0.477 


10.3 


98.89 




53.3d 






12Mg25 


nr, p) 


Na24 


1.368 


99.0 


10.00 




15. Oh 


2.754 


99.0 


12Mg25 


( y, np) 


Na24 


1.368 


99.0 


11.01 




15. Oh 


2.754 


99.0 


19K39 


(V,n) 


K38 


2.170 


100.0 


93.08 




7 . 7mm 






20Ca40 


( y, np) 


K38 


2.170 


100.0 


96.97 




7 . 7mm 






20Ca44 


(y,p) 


K43 


0.374 


84.7 


2.06 




22. 4h 


0.619 


82.4 








1.02 


2.0 


20Ca48 


nr, n) 


Ca47 


1.297 


75.0 


0.18 




4 . 5d 






Ca47 


beta minus 


Sc47 


0.159 


70.0 




decay to 


3 . 4d 
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An unidentified photopeak in the charcoal spectra is 
a 0.706 Mev gamma ray. This peak also appears in the oil 
samples on a regular basis, so it may somehow be associated 
with the presence of the carbon in the samples. However, no 
production reaction to account for the peak could be found. 
A rough half-life is 21 minutes which agrees with the 20.5 
minute half-life of Cll produced in the C12 (Jf,n) Cll 
reaction. 

4. Potting Soil 

A sample of commercial potting soil was analyzed. 
The interst in this sample was the continued development of 
techniques which could be used to analyze petroleum samples. 

Some information on the composition of the soil 
matrix was found by reading the label on the package of soil 
purchased. The contents of the soil included a mixture of 
sandy loam, forest products, charcoal, perlite, iron and peat 
moss. The active ingredients were .05% organic nitrogen and 
.005% thiamin hydrochloride. 

An 8.2 gram sample of soil was irradiated for 10 minutes 
using a 70 Mev electron beam and a 5.3mm lead convertor. The 
radionuclides observed and their probable mechanism of 
production are listed in Table VII. The major and minor 
elements observed were carbon, calcium, magnesium, silicon 
and iron. Unidentified peaks at 0.703 and 0.913 Mev were 
also noted. 
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Table VII. Reactions Observed in the Photoactivation 
of Potting Soil 



Target 

Nucleus 

(abundance) 


Photo- 

reaction 


Residual 
Nucleus 
1/2 Life 


Principal 
Gamma Rays 

Energy 

(Mev) Percent 


14Si29 


( Y* P) 


A128 


1.78 


100.0 


4.70 




2 . 3m 






14Si30 


(}f,P) 


A129 


1.28 


93.0 


3.09 




6 . 6m 






6C12 


( ^/n ) 


Be 7 


0.477 


10.3 


98.89 




53.3d 






12Mg25 


u,p) 


Na24 


1.368 


99.0 


10.00 




15. Oh 


2.754 


99.0 


12Mg26 


( , np 


Na24 


1.368 


99.0 


11.01 




15. Oh 


2.754 


99.0 


19K39 


(iT,n) 


K38 


2.170 


100.0 


93.08 




7.7m 






20Ca40 


U,np) 


K38 


2.170 


100.0 


96.97 


7.7m 






20Ca44 


U,p) 


K43 


0.374 


84.7 


2.06 


22. 4h 


0.619 


82.4 








1.02 


2.0 


20Ca48 


(hn) 


Ca47 


1.297 


75.0 


0.18 




4 . 5d 






20Ca47 


beta minus 


Si47 


0.159 


70.0 




decay to 


3 . 4d 






26Fe57 


( tf,P) 


Mn56 


0.847 


99.0 


2.19 




2. 6h 
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5. Bismuth Germanate (Bi ^ 630 ^ 2 ’) 

A 99.99% pure crystal of BGO was irradiated for 90 
seconds with a direct beam of 100 Mev electrons. The beam 
intensity was measured as 4.8 x 10^ electrons/second. 

Table VIII lists the possible photoreactions for BGO. 
It was assumed that the radionuclides produced would 
predominately be from (V,n), (Y,p) reactions. However, all 
reactions are possible since the irradiating energy will 
produce bremsstrahlung well above the giant dipole resonance. 

Table VIII. Possible Photoreactions with BGO 



Principal 

Gammas 



Target 

Nucleus 

(abundance) 


Photon 

reaction 


Residual 
Nucleus 
Half Life 


Energy 

(Mev) 


Percent 

% 


Percent 

X 

Abundant 


8016 


n 


015 


0.511 


200.0 


200.0 


99.8 




124S 








32Ge70 


np 


Ga 68 


1.078 


4.0 


0.82 


20.5 




68.3m 








32Ge72 


He 3 


Zn69 


0.439 


100.0 


27.4 


27.4 




13. 8h 










np 


Ga70 


1.040 


0.5 


0.14 






21 . 1 m 








32Ge73 




Zn69m 


0.439 


100.0 


7.80 


7.8 




13. 8h 










P 


Ga72 


0.835 


83.3 


6.50 






14. lh 


2.201 


27.3 


2.13 








0.630 


26.7 


2.08 








0.601 


8.2 


0.64 








1.050 


6.9 


0.54 








2.490 


6.8 


0.53 
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Table VIII contd 



Principal 

Gammas 



Target Residual Percent 



Nucleus 

(abundance) 


Photon 

reaction 


Nucleus Energy 

Half Life (Mev) 


Percent 

% 


X 

Abundance 


32Ge74 


o< 


Zn69m 


0.439 


100.0 


36.5 


36.5 




13. 8h 










np 


Ga72 


0.835 


83.3 


30.4 






14. lh 


2.201 


27.3 


9.96 








0.630 


26.7 


9.74 








0.601 


8.2 


2.99 








1.050 


6.9 


2.52 








2.490 


6.8 


2.48 


32Ge74 


P 


Ga73 


0.054 


100.0 


36.5 


36.5 




4 . 9h 


0.295 


99.0 


36.2 








0.74 


6.0 


. 2.19 


32Ge76 


T 


Ga73 


0.054 


100.0 


7.80 


7.8 




4 . 9h 


0.295 


99.0 


7.72 








0.74 


6.0 


0.47 




n 


Ge75 


0.2646 


10.6 


0.83 






82m 








32Ge70 


T 


Ga67 


0.094 


72.12 


14.78 


20.5 




78. 2h 


0.296 


23.03 


4.72 








0.185 


21.95 


4.50 








0.388 


5.74 


1.18 








0.091 


2.20 


0.45 








0.206 


1.52 


0.31 




n 


Ge69 


1.107 


28.10 


5.76 






39. 2h 


0.573 


12.62 


2.59 








0.872 


9.41 


1.93 








1.335 


3.05 


0.63 








0.320 


1.23 


0.25 








0.553 


0.50 


0.10 








0.788 


0.35 


0.07 








0.237 


0.32 


0.07 








1.052 


0.30 


0.06 








1.206 


0.26 


0.05 
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Table VIII. contd. 



Principal 

Ganunas 

Target Residual Percent 

Nucleus Photon Nucleus Energy Percent x 

(abundance) reaction Half Life (Mev) % Abundance 



0.803 99.0 99.0 
0.880 72.0 72.0 
0.516 46.0 46.0 
1.720 36.0 36.0 
0.538 34.0 34.0 
0.343 26.0 26.0 
0.184 21.0 21.0 
0.895 19.0 19.0 
1.099 13.0 13.0 
0.398 10.0 10.0 
1.019 8.0 8.0 
1.596 8.0 8.0 



Table IX is a summary of the photoreactions from 
Table VIII. The information is arranged to assist in the 
identification of the origin of the photopeaks observed in 
the gamma ray spectrum. The last column is five times the 
half-life of the residual nucleus. This is the time when the 
radaioactive nucleus no longer contributes significantly to 
the gamma ray spectrum. 

The irradiated crystal was counted at regular 
intervals for a period of eleven days following the 
irradiation. Figure 9 is a sample spectrum recorded five 
hours after the irradiation. 



83Bi 209 3n Bi206 

100.0 6. 24d 
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Table IX. Summary of Possible Photoreactions of BGO 



Residual 

Nucleus 


Half-Life 


Produced From 


5 x Half-Life 


015 


124S 


016 


10.3m 


Ga 7 0 


21.1m 


Ge72 


105.5m 


Ga68 


68 . 3m 


Ge70 


5 . 7h 


Ge75 


82m 


Ge76 


6 . 8h 


Ga73 


4 . 9h 


Ge74, Ge76 


24. 5h 


Zn69m 


13. 8h 


Ge72, Ge73 , Ge74 


69h 


Ga 7 2 


14. lh 


Ge73, Ge74 


70. 5h 


Ge69 


39. 2h 


Ge70 


196h 


Ga67 


78. 2h 


Ge70 


391h 


3i206 


6 . 24d 


Bi209 


3 Id 



Table X. 


Reactions Observed in the 


Photoactivation of BGO 


Isotope 


Production Reaction 


Method of Verification 


Ga73 


Ge74 (JTfP) Ga73 


.295 photopeak 


Ga 7 2 


Ge73 (3Tfp) Ga72 


Photopeaks plus half- 




and/or 


life of 2.49 Mev 




Ge74 (2f,np) Ga72 


gamma 


Ge69 


Ge70 (Y,n) Ge69 


Photopeaks 


Bi206 


Bi209 (X3n) Bi206 


Photopeaks plus half- 






life of 1.72 Mev gamma 
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A reaction not predicted as probable, but nonetheless 
observed, was the (^3n) reaction of Bi209. This unexpected 
reaction is a good example of the complexity which can be 
introduced into an analysis when higher irradiating energies 
are used. 

There is another interesting way in which the data 
obtained can be manipulated in order to calculate a nuclear 
parameter of interest. That is, the cross section for the 
(tf 3n) reaction of Bi209 can be estimated in the following 
manner . 

First, the photon flux is estimated using equation 
(2) from page 48 and the activity in the 2.49 Mev photopeak 
due to the Ge73 (Y,p) Ga72 rection. Next the activity in the 
1.72 Mev peak due to Bi206 is measured. Finally, the cross 
section for the (^3n) reaction can be calculated by equation 
(2), page 48. This technique of finding the 3i209 (Y/3n) 
cross section will provide a ball park estimate of the cross 
section. For more precise measurements, monochromatic 
photons must be used. Such experiments were done by Berman 
and Fultz [Ref. 25], 

The investigation of BGO concludes the preliminary 
experiments of this thesis. The techniques developed in the 
analysis of charcoal, charcoal residue, potting soil, and BGO 
will be used in the following sections to investigate 
petroleum samples. 
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D. IPAA - PETROLEUM SAMPLES 



1. Background 

The determination of trace elements in oil is a major 
area of interest for the petroleum industry as well as for 
the U.S. Coast Guard. In the former case, the driving force 
is mainly economics. The performance of oil driven 
macahinery, for example, is very sensitive to certain 
concentrations of trace elements such as vanadium. On the 
other hand, the U.S. Coast Guard is interested in trace 
element analysis for its forensic possibilities. 

Various analytical techniques have been used to 
fingerprint oils. These include: IR spectrometry, thin 
layer chromatography, gas chromatography, mass spectrometry, 
atomic absorption, non-disper sive x-ray f lourescence , neutron 
activation analysis, and others. To the best of this 
author's knowledge, no attempts have been made to use photon 
activation analysis to investigate oils. The possibility of 
using instrumental photon activation analysis as a rapid, 
routine method to obtain information about the source of an 
oil spill was the major impetus for these experiments. 

2. Composition of Oil 

In photon activation, it is often helpful to have 
some knowledge of the constituents of a sample which is to be 
irradiated. Such knowledge would be useful in predicting 
which nuclear reactions may occur during an irradiation. It 



65 



would also help in evaluating the experimental data. 
Fortunately, petroleum is a much studied matrix. 



Oil is a generic term which includes crude oil, 
gasoline, lubricating oil, naptha, kerosene and fuel oils. 
Fuel oils have been further classified according to their 
physical properties into five grades: Nos. 1, 2, 4, 5 and 6. 

However, no matter what type of oil is studied, all oils are 
composed primarily of two elements: carbon and hydrogen. 

Sulfur is also present in varying amounts in all crude oils 
and petroleum products. Other minor constituents include 
oxygen and nitrogen compounds. Table XI lists some typical 
concentrations of the major and minor constituents of various 
oils. A general classification for concentration ranges is 
given below. 

1.0 to 100.0% 

0 .< 

io- 



major element 
minor element 
trace element 
subtrace element 



0.01 to 1.0% 

to 10 “ 2 % (i to 100 ppm) 
< 10~^% (< 1 ppm) 



Table XI. Major and Minor Constituents of Various Oils 
[Ref. 18] 





Carbon 


Hydrogen 


Sulphur 


Nitrogen 


Oxygen 








Percent 






Crude oil 


86.06 


13.88 


. 06 


.00 


.00 


Crude oil 


85.05 


12.30 


1.75 


.70 


.00 


Crude oil 


84.00 


12.70 


.75 


1.70 


1.20 


Gasoline 


84.22 


15.73 


.05 


• • 


• • 


Lubricating oil 


85.12 


14.87 


.01 


• • 


• • 


Residual fuel oil 


85.70 


13.93 


.37 


• • 


• • 


Residual fuel oil 


87.10 


12.43 


.47 


• • 


• • 


Residual fuel oil 


86.40 


12. 38 


1.22 


• • 


• • 
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Trace elements which have been identified in oil 



include: Na , Mg, A1 , Si/ P / S/ 01/ K / 0 a / Ti/ V , Or / Mn / Fe / 
Co / Ni / C u / Z n / Q a / As / 3 r / Sr / Mo / Ag , X n / S n / X/ 3a/ Lia/ 
Dy f Au, Pb and Bi. Some of these trace elements are more 
common than others. For example/ Guinn and co-workers 
analyzed 272 samples of oil using neutron activation analysis 
[Ref. 17]. While most of the above elements were found in 
one sample or another, 16 elements were observed on a regular 
basis in most of the samples. Table XII lists these elements 
and their frequency of occurence. 



Table XII. Trace Elements Found Most Frequently in Oil 
[Ref. 17: p. 32] 



Element 


Frequency of Occurence - % 


Na 


99.6 


Br 


99.3 


Cl 


99.3 


V 


91.2 


S 


90.8 


Mn 


89.2 


A1 


83.7 


I 


58.5 


As 


45.2 


Ba 


44.1 


Ni 


41.1 


Gel 


34.1 


Dy 


34.0 


Zn 


33.0 


Co 


18.1 


In 


15.0 



It should be remembered that the trace elements found 
in oils are strongly dependent on the history of the oil. 
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That is to say, the trace elements in oil may not be an 
intrinsic property of the oil. Rather, weathering processes 
such as an oil coming into contact with salt water or sand, 
biological attack, or exposure to the sun may drastically 
alter the trace element composition of the oil. In addition, 
oil inside a holding tank may pick up foreign contaminants 
such as tank scale (rust) or dirt. The implications of the 
above facts of life are that careful attention must be given 
to how samples are collected and even more careful attention 
to how the results of an analysis are interpreted. 

3 . Oil Samples from the Coast Guard Research and 
Development Center 

Five oil samples coded as P-3, P-7, P-11, P-21 and 
P-42 were analyzed. 

a. P-3 Versus P-42 

Samples of P-3 and P-42 were irradiated 
separately for 30 minutes using bremsstrahlung from a 100 Mev 
electron beam and a 5.3mm lead convertor. Each sample was 
counted for 30 minutes at nine hours and at twenty-four hours 
after the irradiation. After a background spectrum was 
subtracted, the peaks observed from both oils were identical. 
At the nine hour mark, an intense peak at .477 Mev, a 
prominent peak at 1.38 Mev and weak peaks at .387, .628 and 
1.94 Mev were observed. The same peaks were present after 
twenty-four hours although with less intensity. The origins 
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of these peaks are listed below. When the spectrua of P-3 
and P-42 were overlapped, there was no significant difference 
between the two spectra. Based on this limited analysis, it 
is concluded that the two oil samples are from the same 
source . 

The peaks observed in samples P-3 and P-24 are 
believed to be from the following reactions. 



Energy (Mev) 


Source 


0.477 


C12 ( ,«*n) Be7 


0.3871 


Ca44 ( y,p)K43 


0.628) 




1.38 


Mg 25 ( 2f,p) Na24 


1.94 


Sum peak 0.477 + 



b. P-7 Versus P-21 

Samples P-7 and P-21 were irradiated separately 
for 30 minutes each using bremsstr ahlung from a 62 Mev 
electron beam and a 5.3mm lead convertor. The average flux 

Q 

for both samples was on the order of 3.5 x 10-° 
photons/cm /sec . 

By increasing the distance of the source to the 
detector and by using a one-half inch lucite disk to 
attenuate positrons, it was possible to count both samples 
one hour after the irradiation. In a fifteen minute count 
the following peaks were observed in the P-21 spectrum. 
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Peak Energy (Mev) Source 



0.511 


Annihilation peak 


1.175") 


C135(Jf,n)Cl34mm 


2.177 b 




3.376J 




1.47 


K40 (background) 


1.67 


sumpeak .511 + 1.175 



No peaks (other than the .511 Mev peak) were observable in 
the P-7 spectrum one hour after the irradiation. 

Additional gamma ray spectrua for both samples 
were recorded at 2, 3, 4, 5 and 10 hours after irradiation. 
Other peaks appearing in the P-21 spectrums were: 



Peak Energy (Mev) 



Source 



0.143 


C135 ( y,n)C134mm 


0.3851 


Ca440f,p)K43 


0.619J 




1.38 l 


Mg25 ( }f,p) Na24 


2.75 ) 




0.708 


Unknown 



Peaks observed in the P-7 spectra were: 



Peak Energy (Mev) 

0.184 

0.511 

0.710 

1.38 



Source 

Unknown 

Annihilation Peak 
Unknown 

Mg25(2T/P) Na24 
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Comparing P-7 and P-21 it is concluded that they 
are from different sources. This observation is based on the 
presence of chlorine and calcium in P-21 but not in P-7, 
c. Irradiation of P-11 

Sample P-11 was irradiated for 30 minutes using 
bremsstrahlung from a 70 Mev electron beam and a 5.3mm lead 
convertor. The flux was estimated as 3 x 10® photons/cm ^/sec. 
The following peaks were observed. 

Peak Energy (Mev) Source 



0.1411 
2.11 ) 


C135(y,n) C134mm 


0.511 


Annihilation peak 


1.02 


Sum peak .511 + .511 


0.703 


Unknown 



This sample appears to be a different oil than 
either P-3, P-7, P-21 or P-42. Because of the presence of 
chlorine, P-11 is different from P-3, P-7 and P-42. Because 
of the absence of any potassium43 peaks, P-11 is different 
from P-21. These results are displayed below. 



- — ____ Oil Sample 
Element ~ — 


P-3 


P-7 


P-11 


P-21 


P-42 


Ca 


X 






X 


X 


Mg 


X 


X 




X 


X 


Cl 






X 


X 





Conclusions: P-3 = P-42 

P-7 4 P-11 4 P-21 4 P-3 or P-42 
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d. Discussion of Results 



After the samples were analyzed, information was 
obtained on the history of the oils. It was learned that 
samples P-3 and P-42 were indeed identical. Furthermore, 
samples P-7 and P-11 were from the same drilling platform but 
were collected one month apart. Finally, sample P-21 was 
from the same general area as samples P-7 and P-11 but was 
obtained from a different platform. 

Considering this a priori knowledge, the 
experimental results, showing samples P-7, P-11, and P-21 are 
from different sources, are reasonable. Although these 
results were based only on an analysis of the minor 
constituents in the oils, the results demonstrate the 
feasibility of using photon activation analysis to 
fingerprint oils. 

4. Oil Samples from the Coast Guard Central Oil Identi- 
fication Laboratory 

Seven samples consisting of 2 crude oils, 2 lube 
oils, 1 fuel oil #5, 1 fuel oil #6 and 1 bunker "C" fuel oil 
were investigated. The irradiation time varied between 10 
minutes and 30 minutes. The electron beam energy was either 
70 or 100 Mev and the lead convertor was used for all 
irradiations . 

The results of these irradiations are summed up with 
one word — uneventful. Except for one lube oil sample where 
the concentration of calcium was determined, the only other 
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significant activity was due to the presence of carbon in the 
oils. These results are not too surprising considering the 
composition of oil as explained previously and considering 
the average beam current for the irradiations was less than 
1 uA. 

Table V lists trace element concentration ranges for 
a sample of fuel oil #6. Although the listed concentrations 
are not necessarily representative of the trace element 
content of fuel oils from different sources, the table does 
imply that calcium may be present in a fuel oil at the minor 
element level. This was indeed the case for one of the lube 
oil samples. 

It was mentioned that the history of an oil may 
drastically alter the concentration of elements found in 
oils. Therefore, it is not uncommon to find certain elements 
in minor amounts in oils. A good example would be a high 
concentration of sodium and chlorine found in an oil spill 
sample collected from the open ocean. 

5 . Automotive, Used Engine Oil 

A sample of used engine oil was irradiated for 10 
minutes using a 100 Mev electron beam and a 5.3mm lead 
convertor. Gamma ray spectra were collected over an eleven 
day period. The gamma rays observed and the probable 
nuclides producing the gamma ray are listed in Table XIII. 
Overlapping gamma rays made it difficult to measure a half- 
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life for most of the photopeaks. The half-life for the 0.278 
Mev gamma ray was measured as approximately 34 hours. 
Analysis of the 0.364 Mev photopeak was complicated by an 
uneven compton distribution from the 0.477 Mev photopeak. A 
rough calculation yielded a half-life of about 60 hours for 
the 0.364 Mev photopeak. A similar observation of an unknown 
gamma ray with an energy of 0.362 Mev and a half-life of 60 
hours was reported by Aras and co-workers in their analysis 
of atmospheric particulate material [Ref. 15: p. 1485]. 



Table XIII. Gamma Rays Observed in the Photoactivation of 
a Sample of Used Engine Oil 



Gamma Ray Energy (Mev) 


Probable Gamma 


Source 


Target Nucleus 


.036 


Unknown 






.078 


Pb X-ray 






.155 


Hgl99 




Pb204 


.185 


Unknown 






.278 


Bal33m, Bal35m 




Bal34 , 135 , 136 , 
137 


.364 


Unknown 






.477 


Be 7 




C12 


.511 


Annihilation 






.620 


Unknown 






.833 


Csl35m 




Bal37 , 138 


.920 


Unknown 






.960 


Pb 2 0 2m 




Pb204 


1.37 


Na24 




Mg25 , 26 


2.74 


Na24 




Mg25 , 26 


6. Determination 


of Calcium in an 


Oil 





Using the methods explained in the theory sections. 



the concentration of calcium in a sample of marine lube oil 
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was determined. In particular/ a comparative assay was 
performed using a lube oil sample and an oil standard spiked 
at 5000 ppm calcium. 

Both the sample and standard were irradiated 
simultaneously for one-half hour using a 70 Mev electron beam 
and a 5.3mm lead convertor. The irradiation and counting 
geometries for the sample and standard were identical. The 
0.619 Mev photopeak from the Ca44(y,p)K43 reaction was 
analyzed because there were no higher energy photopeaks 
appearing in the spectrum whose compton distributions would 
contribute significantly to the observed activity in the 
0.619 Mev photopeak. 

Spectra used for the calculations were obtained at 
around nine hours and twenty-five hours after the 
irradiation. Equation 7 from Chapter II-D is the appropriate 
equation to use. The results follow. 

(Equa. 7) = A :c (tx) e^ t x 

W s Ag(ts)e / ' Ls 

Pertinent Data : 

W g = 1.70 x 10 "^ grams 
A = .0309 hr" 1 
s = 3.4467 grams 



Run 


A x (cpm) 


t x (hr) 


A s (cpm) 


t s (hr) 


1 


27.6 


9.0 


453.1 


9.5 


2 


18.6 


25.0 


264.2 


26.0 
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Calculations: 



Pnn 1 • w = ?7 0309 (9.0) i 7n v i n**2 

W x 453 ,i e *0309(9. b) ^ * ' U X iU 

W x = 1.02 x 10" 3 grams 
C x = W x x 10 6 = 296 ppm 



or C x = W x x 1q2 = *029% 



Run 2 



W x = .18.fi p-° 3 Q 9(25.01^ i.-yo x io-2 

264.2 e .u309uo.0) 



W x = 1.16 x 10 



-3 



grams 



c x = w x/ s x 106 = 336 PP m 
or C x = W x /S x 10 3 = .033% 

Error Analysis : 

The standard deviation of the counting rate is ± Jc~ /t 
where C is the number of counts in a time interval t. It 
is assumed that the half-life of K43 is accurately known as 
22.4 hours. The times of measurements were with ± 1 
minute. Compared to the length of times when the spectra 
were recorded (9, 9.5, 25 and 26 hours), an insignificant 
error is introduced by a ±.1 minute difference. The weights 
of the sample and standard were precisely measured to within 
±..0001 gram. This precision also contributes an 
insignificant error to the analysis. Therefore the largest 
contribution to the error is in the measurement of the 
radioactivity. 
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For run 1 the counting time was 30 minutes. For run 2, 
it was one hour. The counting rates and their standard 
deviation are then. 



Run 


A x (cpm) 


A s (cpm) 


1 


27.6 ± 1.0 


453.1 ± 3.9 


2 


18.6 ± 0.6 


264.2 + 2.1 



The propagation of errors of values that are multiplied 

or divided is related to fractional errors of the 

respective values. In general, for the relationship 

(6-1) W = XY 
Z 



the standard deviation of W is 



( 6 - 2 ) 



$~W = + f 6~x 

W {_ X 



fry 2 
Y 



+ 2 1/2 



In particular, for equation 7, the standard deviations of 



W x is 



(6-3) ( t x ) e W s 

A s (tg)e^ 



M ! *e 



r CHA- )\ 2 



s 



1/2 



When the appropriate numbers are inserted into equation 
(6-3), the final results for the concentration of calcium 
in the sample oil are: 

Run 1 Concentration of Calcium = 296 ± 37 ppm 
Run 2 Concentration of Calcium = 336 ± 38 ppm 

The average of these two values is 316 ± 26 ppm. The 

overall error is then about ± 8% for the measurement. 
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E. MISCELLANEOUS EXPERIMENTS 



1. Spatial Variance of the Photon Flux 

An experiment was conducted in order to obtain some 
information on the spatial distribution of the photon flux. 
Figure 10 shows the experimental set up. 

Copper Wire 



100 Mev 
Electron Beam 

Figure 10. Top View of Experimental Set Up for Flux 
Experiment 

The two targets were glass vials one-half full of 
oil. Each target was wrapped in a sheet of aluminum foil. A 
thin copper wire was placed between the two samples. The 
transverse variance (i.e., x-direction) of the flux was first 
investigated by doing a gross count of aluminum foils A and 
B. Using the procedures as explained in Section II-C-1, the 
average value of the flux for foil A was 2.27 x 10® 
photons/cm Vsec. For foil B the flux was 0.27 x 10® 
photons/cm ^/sec. This large difference in average flux is 
attributed to misalignment of the targets. 

The Z-dependence of the photon flux was determined by 
cutting foil A into six sections. Each section was about 1 



7 



B 



.6.35 cm. 
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cm long. The sections were then counted individually and the 
flux determined. The results were: 



Section 

1 

2 

3 

4 

5 

6 



Flux (photons/cm^/sec) 



6.58 x 10 
3.33 x 10 



3.66 x 10? 



2.73 x 10 
2.41 x 10 
3.31 x 10 



When the results are normalized to the flux for section 1, 
the Z-dependence of the flux is given by the following 
histogram. 



4 > 



relative 

units 




Distance in Z-direction (cm ) 



One explanation for the above profile is that the target was 
completely wrapped with aluminum foil. Therefore, the end 
sections expose more surface area to the photon flux than the 
foil on the sides. 

As a check on the observed flux profile above, a 
similar calculation was done using the copper wire. The 
reaction used was the Cu65 ( ^,n)Cu64 reaction. The half-life 
of Cu64 is 12.25 hours and it decays by positron emission and 
electron capture. Thirty-eight percent of the time a 0.511 
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gamma ray is emitted. The 0.511 gamma ray was used for 
calculations. The results were: 



Section Flux 

1 

2 

3 

4 

5 



(photons/cm /sec) 

3.54 x 10 9 

2.21 x 10 9 

1.21 x 10 9 
.60 x 10 9 
.32 x 10 9 



When these results are normalized, the flux profile is given 



by: 



relative 

units 




This is a different flux profile than that calculated 
using the aluminum foil. This is not surprising since only 
relative calculations were done rather than absolute calcula- 
tions. For example, no nuclear data for the (Jf,He^) cross 
section for the aluminum calculation was available. A value 
of 1 millibarn was used. Assuming accurate values are known 
for all variables in the calculation (i.e., cross sections. 
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detector efficiencies, decay schemes, etc.)/ it would be 
possible to do an absolute determination of the flux. 



The bottom line of this experiment is that it is very 
difficult to reproduce the same photon flux in different runs 
of an experiment. Therefore, as mentioned previously, the 
comparative technique is most often used in activation 
analysis. 

2 . IPAA of Oils Spiked With Calcium and Iron 

A sample of oil containing 5000 ppm iron and a sample 
of oil containing 5000 ppm calcium were irradiated. The 
purpose of the irradiation was to observe the gamma rays 
emitted due to a known amount of element in a sample. The 
spectra obtained would then be useful in identifying gamma 
rays due to either iron or calcium in an unknown sample. 

Table XIV lists the reactions observed in the calcium 
sample. Table XV lists the reactions observed in the iron 
sample. In both samples a strong peak at 0.710 Mev was 
observed. In neither case could this peak be attributed to 
the iron or calcium activity. 
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Table XIV. Reactions Observed in the Photoactivation of 
an Oil Sample Spiked with 5000 ppm Calcium 



Principal 
Gamma Rays 

Target Residual 

Nucleus Photo- Nucleus Energy 

(abundance) reaction 1/2 Life (Mev) Percent 



20Ca43 


(r,2p) 


0.145 


or,p) 


20Ca44 

2.06 


( y f He3) 




( y,np) 




<y,P) 


20Ca48 

0.18 


U, n) 


20Ca47 


beta minus 
decay to 



M > 
• h 
00 ** 
tr h 


1.293 


99.2 


K42 
12. 4h 


1.528 


18.0 


Ar 41 
1 . 8h 


1.293 


99.2 


K42 
12. 4h 


1.528 


18.0 


K43 


0.374 


84.7 


22. 4h 


0.619 


82.4 




1.02 


2.0 


Ca47 
4 . 5d 


1.297 


75.0 


Sc47 
3 . 4d 


0.159 


70.0 
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Table XV. Reactions Observed in the Photoactivation of an 
Oil Sample Spiked with 5000 ppm Iron 









Principal 








Gamma 


Rays 


Target 




Residual 






Nucleus 


Photo- 


Nucleus 


Energy 




(abundance) 


reaction 


1/2 Life 


(Mev) 


Percent 


26Fe54 


(Y,2n) 


Fe52 


0.165 


100.0 


5.82 




8 . 2h 






26Fe56 


(ar, np) 


Mn54 


0.834 


100.0 


91.66 




312. 5d 






26Fe57 


(y»p) 


Mn56 


0.847 


99.0 


2.19 




2 . 6h 


1.811 


29.0 








2.110 


15.5 


26Fe58 


( 3r",np) 


Mn56 


0.847 


99.0 


0.33 




2 . 6h 


1.811 


29.0 








2.110 


15.5 
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IV. SUMMARY AND CONCLUSIONS 



The first part of this thesis reviewed the basic 
principles of photon activation analysis. It was noted that 
the technique had been used successfully in the investigation 
of a variety of materials. In the second part of the thesis, 
experiments with general samples showed that the major and 
minor elements in a sample could be easily activated. The 
most common elements observed in biological and environmental 
samples were carbon, silicon, calcium, potassium, magnesium 
and iron. 

After investigating the general samples, various samples 
of oils were analyzed. Even though no trace elements were 
observed, it was still possible to make tentative comparisons 
of oils based on the minor constituents. A quantitative 
analysis for calcium showed that it was possible to detect 
elements at a concentration of .03% (300 ppm). 

An interesting peak appeared in many spectra at 0.710 
Mev. It is believed that this peak is a sum peak due to the 
intense activity of the 0.511 Mev peak. For example, if the 
0.511 Mev peak and its backscatter peak are counted 
simultaneously, a peak would appear around 0.7 Mev. The 
half-life associated with the 0.710 Mev peak, measured as 21 
minutes, is further evidence for this hypothesis. It is known 
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that the 0.511 peak is mainly due to carbon 11 which has a 
half-life of 20.5 minutes. 

A number of conclusions can be stated based on the 
experimental results. First, even with the low beam current 
of the LINAC, it is possible to determine major and minor 
elements in a matrix. A lower limit of detection for most 
elements is estimated as 100 ppm. Since the range of trace 
elements is 1 to 100 ppm, higher beam currents are required 
for trace element analysis. In addition, if a complex matrix 
such as oil is to be analyzed, many gamma rays can be 
expected when the sample is counted. Therefore, to avoid 
overlapping gamma rays, high resolution is important and 
hence a Ge{Li) detector would be required for any 
quantitative results. 

A second conclusion of this thesis is that it would be 
beneficial to continue to investigate oil using photon 
activation analysis. One reason to continue is purely 
academic. That is, photon activation is a highly sensitive 
technique which has not been used previously to fingerprint 
oil spills. The possibility of identifying new elements not 
previously reported in oils is in itself sufficient reason 
for further analysis. Perhaps a more practical reason is a 
recent oil spill case which occurred in the 13th Coast Guard 
District. In that incident, samples were collected from the 
spill and suspect vessels and a positive match was made by 
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the Coast Guard's Central Oil Identification Laboratory. The 
methods used included infra-red spectrometry, gas 
chromatography, thin layer chromatography, and uv 
f lourescence. Owners of the suspect vessel also commissioned 
an analysis of the oil samples using neutron activation 
analysis. The interpretation of the results of the neutron 
activation analysis apparently did not agree with the Coast 
Guard's tests. Armed with this knowledge, there is the 
possibility that the suspect vessel's owners may challenge 
the validity of the Coast Guard tests. Even though the Coast 
Guard tests have been upheld by the courts, the fact that 
activation analysis is being considered as a technique to 
perhaps discredit the present tests is reason for concern. 
To anticipate any problems in this area, it is recommended 
that the Coast Guard consider activation analysis as an 
additional technique of analyzing oils. For example, in oil 
spill cases where large amounts of money are involved for 
cleanup and fines, it may be wise to contract for an 
activation analysis. Towards this end, it would also be 
helpful to have more than one Coast Guard officer familiar 
with activation analysis procedures. One way to do this 
would be to collaborate with the LINAC branch of the National 
Bureau of Standards. One or two weeks at NBS would be 
sufficient time to generate enough data for weeks of 
analysis. The study could easily be developed into an 
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which could involve 



interdisciplinary thesis project, 
students from the physics, computer science and operations 
research curricula. What would be needed to initiate such a 
project would be the sponsorship of an appropriate command 
such as COIL. 

As the first step of any possible further work with oil 
and photon activation analysis. Appendix A is a theoretical 
analysis of oil using a 125uA beam current and a 60 cm^ 
Ge(Li) detector. 

The final conclusion of this thesis is that the analysis 
of oil for trace elements using the Naval Postgraduate 
School's LINAC is not feasible. However, elements in 
concentrations down to .01% (100 ppm) may be detected. One 
possible way to observe elements in oil other than carbon is 
to burn the oil prior to irradiation. This will eliminate 
the carbon interference and concentrate the trace elements in 
the remaining ash. This ash almost always contains amounts of 
aluminum, calcium, iron, magnesium, nickle, sodium, silicon 
and vanadium ranging from 0.015% to 0.05% of the ash 
(Ref. 18: p. 94]. The problem with this method is that the 
maximum amount of ash in oil is about 0.2%. It would 
therefore take 500 grams of oil to produce a 1 gram sample of 
ash. 
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APPENDIX A 



A THEORETICAL ANALYSIS OF A FUEL OIL USING PHOTON ACTIVATION 

ANALYSIS 

Trace element analysis is not possible with the NPGS 
linear accelerator. Therefore, a hypothetical experiment 
will be described which will use photon activation analysis 
to analyze an oil sample. The sample to be studied is a 
number 6 residual fuel oil. The concentration of 26 trace 
elements in this oil were determined in an inter laboratory 
comparison study [Ref. 19: p. 244]. These elements and their 
concentrations are listed in Table A-l. The concentration 
values underlined were determined by neutron activation 
analysis . 

The experiment to be discussed can be considered a "what 
if" experiment. That is to say, "What if it is assumed that 
the known trace element concentrations in an oil sample are 
in fact unknown. What trace elements could then be detected 
in a photon activation analysis?" A convenient analysis for 
this purpose would be one similar to the experiment done by 
Chattopadhyay and Jervis in their photon activation analysis 
of market-garden soils [Ref. 12]. They irradiated a one gram 
sample of soil for a standard irradiation time of one hour. 
The electron beam energy varied from 8 to 44 Mev. A 3 mm 
thick, water-cooled tungsten convertor provided an average 
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flux of 2 x 1C)13 photons/cm ^/sec. The sample was placed 1 cm 
away from the rear face of the convertor. The typical 
average beam current measured on the convertor ranged from 35 
to 190 uA. Samples were allowed to "cool" for one hour after 
the irradiation, then the gamma rays were observed with a 60 
cm^ Ge(Li) detector coupled with a 4096 channel multi-channel 
analyzer. The FWHM of the detector was 1.9 KeV at the 1.332 
Mev photopeak of Co60. The efficiency of the detector was 
measured as 6.7% with respect to a standard Nal(Tl) detector. 

This analysis of the fuel oil sample will assume all the 
conditions above except only one irradiation energy will be 
used. The energy chosen is 25 Mev implying a beam current of 
125 uA. The reason for this choice is that Chattopadhyay and 
Jervis feel that a lower electron energy and a higher beam 
current are sufficient to estimate a number of elements. 
Using a higher irradiation energy will produce more activity, 
but it will also enhance interfering reactions. 

Table A-2 lists the probable photoreactions expected for 
the elements in Table A-l. The considerations in deciding 
which reactions are probable are: (a) The (V,n), ( V/P) 
reaction cross sections are about equal up to atomic number 
30. Above Z =30 the (}f,n) reaction will predominate, 
(b) Since the irradiation energy is 25 Mev, reactions having 
a threshhold energy above 20 Mev are not likely to produce 
any measurable activity. (c) The irradiation time of one 
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hour will favor nuclear products with half-lives less than 20 
minutes. (d) However, by waiting one hour before counting, 
any radionuclide with a half-life less than 10 minutes will 
probably not be detected. 

Table A- 3 lists the photoreactions of Table A-2 in order 
of increasing gamma ray energy. Interference-free gamma rays 
have been identified with an asterisk. Gamma rays falling 
within ±4 KeV of each other were considered to be mutually 
interfering under the counting system used in this analysis. 

The interference-free gamma rays in Table A-3 can be used 
to identify an element in the sample. This is done by 
measuring a half-life corresponding to each photopeak. Using 
the gamma ray energy and half-life information, a 
radionuclide can be associated with each photopeak. The 
target nucleus and therefore an element in the oil sample can 
then be deduced. A summary of the elements which can be 
determined in this manner is listed below. The number in 
parenthesis is the number of interference-free gamma rays 
associated with each element. 



Elements Which Can be Detected in Oil 



(Based on interference-free gammas only) 



Ag (4) 
Zn (3) 
Cd (3) 
Fe(3) 



As (2) 
Sb ( 2) 
Se(2) 
Cu ( 2) 



V (2) 
Mg (2) 
Mn (1) 
Hg (1) 



Ca(l) 
Ni (1) 
Cr (1) 
Pb (1) 
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For some of the elements above which yield only one 
interference-free gamma ray, it is possible to take advantage 
of different decay rates to obtain additional interference- 
free gammas. For instance, if the sample is counted 3 days 
after the irradiation, the 0.065 Mev gamma of Snll9m will be 
free from any interference from shorter lived isotopes. 
Another example is the determination of lead. The photopeaks 
of Zn63 and Pb204m overlap at 0.961 Mev. However, the half- 
lives of Zn63 and Pb204m are 38 minutes and 3.6 hours 
respectively. Thus, by counting the sample at around 5 hours 
after the irradiation, the activity in the 0.961 Mev peak 
will be due to lead alone. Finally, the 0.374 Mev peak of 
K43 becomes interference-free if one counts the sample after 
the activity of Hgl99m decays to an insignificant level. 

A factor not considered in the above discussion is the 
relative abundance of the trace elements in the oil sample. 
This was because the assumption was made that the trace 
elements in the oil were unknown. However, if the 
concentration of an element is known, the question must be 
asked is it within the detection limits of the analysis. The 
detection limits for this analysis were assumed to be 
identical to the detection limits determined by Chattopahyay 
and Jervis. This assumption is not totally valid because 
detection limits must be restricted to a particular matrix 
studied and the concentrations of interfering elements in 
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that matrix. However, since the experimental conditions 
assumed for the oil sample were identical to the conditions 
used to analyze the soil sample, it is reasonable to use the 
same detection limits for both samples. When these detection 
limits are applied to the oil analysis, it is found that 
manganese, selinium, silver, antimony and mercury cannot be 
detected. The lack of gamma .rays due to these elements does, 
however, allow the presence of barium to be determined. The 
list of elements which can then be determined in the oil 
sample becomes: Zn, Cd, Fe, As, Cu, V, Mg, Sn, Ca, Ni, Cr, 
Pb and Ba. It is interesting to compare the above list to a 
list of the elements determined by neutron activation 
analysis. Elements which can be determined by either method 
include: Ca, V, Cr, Fe, Ni, Zn and As. Elements which can 
be determined by NAA but not PAA include: Sb, Hg, Ag, Se, Mn 
and Na. Finally, elements detectable by PAA but not NAA 
include: Mg, Cu, Cd, Sn, Ba and Pb. In both cases 13 
elements are detectable. However, the use of PAA allows 
different elements to be determined than those determined 
with NAA and vice versa. This comparison, of course, only 
applies to this particular example. Nonetheless, the general 
conclusion that PAA and NAA will produce different nuclides 
is valid. 

Another point to consider in this analysis is the 
presence of trace elements in the oil sample other than the 
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elements listed in Table A-l. For example, Al, Cl, Co, Ga, 
Br, In, I and Dy are also commonly found in oil. (See Table 
XII). It is not known whether or not these elements were 

detected in the oil sample reported in Reference 19. For the 

/ 

purpose of this analysis it was assumed that they were not 
present in the oil. 

The last comment concerning this study is that for the 
same oil sample, a different set of experimental conditions 
will yield different results. Therefore, by varying the 
irradiating energy, the beam current, or the irradiation 
time, certain reactions can be enhanced and other reactions 
can be diminished or eliminated. This flexibility of 
optimizing certain reactions is one of the major advantages 
of photon activation analysis. 
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Table A— 1 . 


Trace Element Concentrations of a Residual 
Fuel Oil Number 6 [Ref. 19: p. 244] 


Elements 


Concentration (ppm) 


Li 

Be 

B 

F 

Na 

Mg 

Si 

K 

Ca 

V 

Cr 

Mn 

Fe 

Ni 

Cu 

Zn 

As 

Se 

Sr 

Ag 

Cd 

Sn 

Sb 

Ba 

Hg 

Pb 


3 

0.5 

0.2 

0.004 

7.7 

3 
30 

5 

400 

113 

1 

1 

20 

62 

1 

1.4 

0.35 

0.15 

0.5 

0.02 

1 

5 

0.01 

5 

0.01 

4 
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Table A-2. Probable Photoreactions for Sample Fuel Oil 



Target 

Nucleus 

(abundance) 



Photo- 

reaction 

(E t ) 



Residual 

Nucleus 



1/2 Life 



Principal 

Gammas 

Intensity 

Energy Intensity times 
(Mev) Percent Abundance 



Cll 0.511 200.0 199.0 

20.5m (annihilation peak) 



6C12 


n 


98.89 


18.7 


12Mg25 


P 


10.0 


12.1 


20Ca44 


P 


2.06 


12.2 


20Ca48 


n 


0.18 


9.9 



23V50 

0.24 


2n 

20.9 


24Cr 50 
4.35 


n 

12.9 


24Cr 52 
83.79 


n 

12.0 


25Mn55 

100.0 


n 

10.2 


26Fe57 

2.19 


P 

10.6 


28Ni58 

67.77 


n 

12.1 


28Ni 62 
3.66 


P 

11.1 



Na24 


1.368 


15. Oh 


2.754 


K43 


0.374 


22. 4h 


0.619 


Ca47 
4 . 5d 


1.297 


B" to 
3 . 4d 


Sc470 . 159 


V48 


0.983 


16. Id 


1.31 


Cr49 


0.091 


41.9m 


0.062 

0.153 


Cr 51 
27. 8d 


0.320 


Mn54 
312. 5d 


0.834 


Mn56 


0.847 


2 . 6h 


1.811 

2.110 


Ni57 


0.127 


36h 


1.378 

1.919 


Co61 

1.7h 


0.067 



99.0 


9.9 


99.0 


9.9 


84.7 


1.8 


82.4 


1.7 


75.0 


0.1 


70.0 


0.1 


100.0 


0.2 


97.0 


0.2 


30.0 


1.3 


15.0 


0.7 


14.0 


0.6 


9.9 


8.3 


100.0 


100.0 


99.0 


2.2 


29.0 


0.6 


15.5 


0.3 


14.0 


9.5 


86.0 


58.2 


14.0 


9.5 


100.0 


3.7 
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Table A- 2 Contd 



Target 

Nucleus 

(abundance) 



Photo- Residual 
reaction Nucleus 
(E t ) 1/2 Life 



Principal 

Gammas 



Energy Intensity 
(Mev) Percent 



29Cu63 


2p 


Co61 


0.067 


100.0 


69.17 


17.2 


1 . 7h 








2n 


Cu61 


0.284 


12.0 




19.7 


3 . 4h 


0.655 


11.0 








0.674 


4.0 








0.38 


2.5 








0.58 


1.5 








0.94 


1.4 


29Cu65 


c* 


Co61 


0.067 


100.0 


30.83 


6.7 


1 . 7h 






30Zn64 


n 


Zn63 


0.669 


8.0 


48.89 


11.9 


38.4m 


0.962 


6.0 


30Zn66 


n 


Zn65 


1.115 


50.6 


27.81 


11.0 


244d 






30Zn68 


P 


Cu67 


0.184 


40.9 


18.57 


10.0 


61. 6d 


0.093 


24.5 


30Zn70 


n 


Zn69m 


0.439 


100.0 


0.62 


9.2 


13. 8h 






33As75 


n 


As 7 4 


0.596 


61.8 


100.0 


10.3 


17. 7d 


0.635 


14.0 


34Se74 


n 


Se73 


0.359 


100.0 


0.87 


12.1 


7 . lh 


0.066 


100.0 


34Se76 


n 


Se75 


0.265 


59.5 


9.02 


11.2 


120. 4d 


0.136 


57.3 








0.279 


25.2 








0.121 


16.9 








0.400 


12.4 


34Se77 


P 


As 7 6 


0.559 


42.0 


7.58 


9.6 


26. 3h 







Intensity 

times 

Abundance 



69.2 



8.3 

7.6 

2.8 

1.7 

1.0 

1.0 

30.8 



3.9 

2.9 

14.2 



7.6 

4.6 

0.6 



61.8 

14.0 

0.9 

0.9 

5.4 

5.2 

2.3 

1.5 

1.1 

3.2 
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Table A-2 Contd 



Target 

Nucleus 

(abundance) 



Photo- 

reaction 



(E t ) 



Residual 

Nucleus 



1/2 Life 



Principal 

Gammas 

Intensity 

Energy Intensity times 
(Mev) Percent Abundance 



34Se82 


n 


Se81m 


9.19 


9.3 


57m 


38Sr 86 


n 


Sr 85m 


9.87 


11.5 


70m 


38Sr 87 


2n 


Sr 85m 


7.04 


19.9 


70m 




V 1 


Sr87m 




M4 


2 . 8h 


38Sr 88 


n 


Sr87m 


82.53 


11.1 


2 . 8h 


47Agl07 


n 


Agl06m 


51.82 


9.5 


8 . 4d 



48Cdll2 


n 


Cdlllm 


24.07 


9.4 


48.6 


48Cdll3 


2n 


Cdlllm 


12.26 


15.9 


48.6 


48Cdll6 


n 


Cdll5 


7.58 


8.7 


53. 5h 


50Snll8 


n 


Snll7m 


24.03 


9.3 


14. Od 



0.10 3 


100.0 


9.2 


0.237 

and 

0.231 


86.0 


8.5 


0.237 

and 

0.231 


86.0 


6.1 


0.388 


99.0 


7.0 


0.388 


99.0 


81.7 



0.511 


89.5 


46.4 


1.046 


31.1 


16.1 


0.808 


25.0 


13.0 


0.717 


23.3 


12.1 


0.406 


22.3 


11.6 


0.748 


22.1 


11.5 


0.616 


20.3 


10.5 


1.528 


20.2 


10.5 


+ 13 


more gammas 




0.247 


100.0 


24.1 


0.150 


100.0 


24.1 


0.247 


100.0 


12.3 


0.150 


100.0 


12.3 


0.335 


95.0 


7.2 


0.53 


26.5 


2.0 


0.159 


100.0 


24.0 


0.158 


100.0 


24.0 
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Table A-2 Contd 



Target 

Nucleus 

(abundance) 



Photo- 

reaction 



(E t ) 



Residual 
Nucleus 
1/2 Life 



Principal 

Gammas 

Intensity 

Energy Intensity times 
(Mev) Percent Abundance 



50Snl20 


n 


32.85 


9.1 


50Snl24 


n 


5.94 


8.5 


51Sbl21 


n 


57.25 


9.3 



5lSbl23 


n 


42.75 


8.9 


56Bal34 


n 


2.42 


9.3 


56Bal35 


V 1 


6.59 


M4 


56Bal36 


n 


7.81 


9.2 


80Hgl98 

10.02 


n 

8.6 


80Hg200 

23.13 


n 

8.0 


80Hg201 

13.22 


2n 

14.3 


82Pb204 


n 


1.48 


8.2 




2n 

15.2 



Snll9m 

245d 


0.065 


Snl23m 

40m 


0.160 


Sbl20 

15.9m 


1.171 


Sbl20 
5 . 8d 


1.170 

0.200 


Sbl22 
2 . 7d 


0.564 

0.686 


Bal33m 
38. 9h 


0.276 


Bal35m 
28. 7h 


0.268 


Bal35m 
28. 7h 


0.268 


Hgl97 
64. lh 


0.077 


Hgl99m 

43m 


0.375 

0.158 


Hgl99m 

43m 


0.375 

0.158 


Pb203 
52. lh 


0.279 


Pb202m 
3 . 6h 


0.961 

0.422 

0.658 

0.129 



100.0 


32.9 


100.0 


5.9 


1.3 


0.7 


100.0 


57.3 


100.0 


57.3 


66.3 


28.3 


3.3 


1.4 


100.0 


2.4 


100.0 


6.6 


100.0 


7.8 


99.9 


10.0 


100.0 


23.1 


100.0 


23.1 


100.0 


13.2 


100.0 


13.2 


99.2 


1.5 


90.0 


1.3 


90.0 


1.3 


78.2 


1.2 


78.2 


1.2 
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Table A-3. Probable Photoreactions for Sample Fuel Oil 
(Listed in Order of Increasing Gamma Ray 
Energy) 



Intensity 

times 

Abundance 



Energy 


Target 


Reaction 


Product 


Half-Life 


% 


0.062 


Cr 50 


n 


Cr49 


41.9m 


0.7 


0.065 


Snl20 


n 


Snll9m 


245d 


32.9 


0.066 


Se74 


n 


Se73 


7 . lh 


0.9 


0.067 


Ni62 


p 


Co61 


1. 7h 


3.7 




Cu63 


2p 


Co61 


1 . 7h 


69.2 




Cu65 




Co61 


1. 7h 


38.8 


0.077* 


Hgl98 


n 


Hgl97 


64. lh 


10.0 


0.091 


Cr50 


n 


Cr49 


41.9m 


11.3 


0.093 


Zn68 


P 


Cu67 


61. 6h 


4.6 


0.103* 


Se82 


n 


Se81m 


57m 


9.2 


0.121 


Se76 


n 


Se75 


120d 


1.5 


0.127 


Ni58 


n 


Ni57 


36h 


9.5 


0.129 


Pb204 


2n 


Pb202m 


3 . 6h 


1.2 


0.136 


Se76 


n 


Se75 


120d 


5.2 


0.150 


Cdll2 


n 


Cdlllm 


48.6m 


24.1 


0.153 


Cr 50 


n 


Cr 49 


41.9m 


0.6 


0.158 


Hg200 


n 


Hgl99m 


43m 


23.1 


0.159 


Snll8 


n 


Snll7m 


14d 


24.0 


0.160 


Snl24 


n 


Snl23m 


40m 


5.9 


0.184* 


Zn6 8 


P 


Cu67 


61. 6h 


7.6 


0.200* 


Sbl21 


n 


Sbl20 


5. 8d 


57.3 


0.231 


Sr 86 


n 


Sr 85m 


70m 


6.1 




Sr 87 


2n 


Sr85m 


70m 


6.1 


0.237 


Sr 86 


n 


Sr 85m 


70m 


6.1 




Sr 87 


2n 


Sr 85m 


70m 


6.1 


0.247* 


Cdll2 


n 


Cdlllm 


48 . 6m 


24.1 
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Table A-3 Contd 



Intensity 

times 

Abundance 



Energy 


Target 


Reaction 


Product 


Half-Life 


% 


0.265 


Cdll3 


2n 


Cdlllm 


48 . 6m 


12.3 




Se76 


n. 


Se75 


120d 


5.4 


0.268 


Bal35 


1 


Bal35m 


28. 7h 


6.7 




Bal36 


n 


Bal35m 


28. 7h 


7.8 


0.276 


Bal34 


n 


Bal35m 


38. 9h 


2.4 


0.279 


Hg204 


n 


Hg203 


46d 


6.7 




Pb204 


n 


Pb203 


52h 


1.5 




Se76 


n 


Se75 


120d 


2.3 


0.284 


Cu63 


2n 


Cu61 


3 . 4h 


8.3 


0.320* 


Cr 52 


n 


C451 


27. 8d 


8.3 


0.335* 


Cdll6 


n 


Cdll5 


53. 5h 


7.2 


0.359* 


Se74 


n 


Se73 


7. Ih 


0.9 


0.374 


Ca44 


P 


K43 


22. 4h 


1.8 


0.375 


Hg200 


n 


Hgl99m 


43m 


23.1 




Hg201 


2n 


Hgl99m 


43m 


13.2 


0.38 


Cu63 


2n 


Cu61 


3. 4h 


1.7 


0.388 


Sr 87 


1 


Sr87m 


2. 8h 


7.0 




Sr 88 


n 


Sr87m 


2 . 8h 


81.7 


0.400 


Se76 


n 


Se75 


120d 


1.1 


0.406 


Agio 7 


n 


AglO 6m 


8. 4d 


11.7 


0.422* 


Pb204 


2n 


Pb 2 0 2m 


3 . 6h 


1.3 


0.439* 


Zn70 


n 


Zn69m 


13. 8h 


0.6 


0.511 


C12 


n 


Cll 


20.5m 


annihila- 
tion peak 


0.511 


AglO 7 


n 


Agl06m 


8 . 4d 


46.4 


0.514 


Sr 86 


n 


Sr 85 


64. 5d 


9.9 


0.53* 


Cdll6 


n 


Cdll5 


53. 5h 


2.0 


0.559 


Se77 


P 


As 7 6 


26. 3h 


3.2 


0.564 


Sbl23 


n 


Sbl22 


2 . 7d 


28.3 
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Table A-3 Contd 



Intensity 

times 

Abundance 



Energy 


Target 


Reaction 


Product 


Half-Life 


% 


0.58* 


Cu63 


2n 


Cu61 


3 . 4h 


1.0 


0.596* 


As 7 5 


n 


As74 


17. 7d 


61.8 


0.616 


Agl07 


n 


AglO 6m 


8 . 4d 


10.5 


0.619 


Ca44 


p 


K43 


22. 4h 


1.7 


0.635* 


As75 


n 


As 7 4 


11 .Id. 


14.0 


0.655 


Cu63 


2n 


Cu61 


3 . 4h 


7.6 


0.658 


Pb204 


2n 


Pb202m 


3 . 6h 


1.2 


0.669 


Zn64 


n 


Zn63 


38m 


3.9 


0.674 


Cu63 


2n 


Cu61 


3 . 4h 


2.8 


0.686* 


Sbl23 


n 


Sbl22 


2. 7d 


1.4 


0.717* 


AglO 7 


n 


AglO 6m 


8 . 4d 


12.1 


0.748* 


Agio 7 


n 


AglO 6m 


8 . 4d 


11.5 


0.808* 


AglO-7 


n 


AglO 6m 


*0 

rr 

• 

00 


13.0 


0.834* 


Mn55 


n 


Mn54 


312d 


100.0 


0.847* 


Fe57 


P 


Mn56 


2 . 6h 


2.2 


0.94* 


Cu63 


2n 


Cu61 


3 . 4h 


1.0 


0.961 


Pb204 


2n 


Pb202m 


3 . 6h 


1.3 


0.962 


Zn64 


n 


Zn63 


38m 


2.9 


0.983* 


V50 


2n 


V48 


16. Id 


0.2 


1.046* 


AglO 7 


n 


AglO 6m 


8 . 4d 


16.1 


1.115* 


Zn66 


n 


Zn65 


243d 


14.2 


1.170 


Sbl21 


n 


Sbl20 


5 . 8d 


57.3 


1.171 


Sbl21 


n 


Sbl20 


15.9m 


0.7 


1.297* 


Ca48 


n 


Ca47 


4 . 5d 


0.1 




B - decay to 


Sc 4 7 

(.159 Mev) 


3 . 4d 


0.1 


1.31* 


V50 


2n 


V48 


16. Id 


0.2 



101 



Table A- 3 Contd 



Energy 


Target 


Reaction 


Product 


Half-Life 


Intensity 

times 

Abundance 

% 


1.369* 


Mg 2 5 


P 


Na24 


15. Oh 


9.9 


1.378* 


Ni58 


n 


Ni57 


36h 


58.2 


1.528* 


Agio 7 


n 


Agio 6m 


8 . 4d 


10.5 


1.811* 


Fe57 


P 


Mn56 


2. 6h 


0.6 


1.919* 


Ni58 


n 


Ni 57 


36h 


9.5 


2.110* 


Fe57 


P 


Mn56 


2. 6h 


0.3 


2.754* 


Mg 25 


P 


Na24 


15. Oh 


9.9 



* Interference free reactions 
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